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Abstract. Here we aim at justifying rigorously different types of physically relevant diffu¬ 
sive limits for radiative flows. For simplicity, we consider the barotropic situation, and adopt 
the so-called PI-approximation of the radiative transfer equation. In the critical functional 
framework, we establish the existence of global-in-time strong solutions corresponding to 
small enough data, and exhibit uniform estimates with respect to the coefficients of the sys¬ 
tem. Combining with standard compactness arguments, this enables us to justify rigorously 
the convergence of the solutions to the expected limit systems. 

Our results hold true in the whole space R" as well as in a periodic box T" with n > 2. 

Keywords: Radiation hydrodynamics, Navier-Stokes system, diffusive limit, critical reg¬ 
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1. Introduction 

We consider the barotropic version of a model of radiation hydrodynamics. Our main goal 
is to provide the rigorous justification of asymptotics that have been investigated formally 
and numerically by Lowrie, Morel and Hittinger m. and mathematically by the second 
author and S. Necasova in uni [mils] in the finite energy weak solutions framework. 

The fluid is described by standard classical fluid mechanics for the mass density g and 
the velocity field u as functions of the time t G M+ and of the (Eulerian) spatial coordinate 
a: € D where D is either the whole space M"" or some periodic box T” with n > 2. 
Radiation acts through some radiative momentum source Sp which is given by 

Sp = - r f uSdudu, 

C Jo Js^-^ 

where c is the light speed. 

The radiative source S = S{t,x,Ld,p) depends on the direction vector uj G 5"’“^ (where 
denotes the unit sphere of M”), and on the frequency u > 0 of the photons, and is 
given by 

S = aa{B{p, g) — + as{i — where I := _ j I du. 

The radiative intensity I obeys the transfer equation 

(1.1) - dtX + u-V:cI= S in (0,r) X D X 5"-^ X (0,cx)). 

c 
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In the present paper, as in mu, we make the following simplifying assumptions 

(1) Isotropy : the transport coefficients da and ag are independent of w; 

(2) ‘Gray’ hypothesis : da and d^ are independent of ly; 

(3) ‘PI hypothesis’ : the averaged radiative intensity I := I dv is given by the 
ansatz 

(1.2) I = Iq u) ■ Ii, 
where Iq and Ii are independent of uj and 

Plugging (jl.2l) in (jl.Ij) . and computing the 0th and 1st order momentum with respect to 
oj, we find out the following evolution equations for Iq and Ii (keeping the same notation 
B for the distribution function averaged in v) 

(1.3) - dflo + - diva;/i = aa{g)iB{g) - Iq), 

c n 


(1-4) - dth + Vxh = {<Ta{g) +aa{g))Ii. 

c 

Besides, the radiative force is now given by 

.5) s^=('^£M±^)r„ 

In order to identify the most relevant asymptotic regimes, we rewrite the equations in 
dimensionless form. To this end, introduce some reference hydrodynamical quantities (length, 
time, velocity, density, pressure): Z, T, U, g, p, and reference radiative quantities (radiative 
intensity, absorption and scattering coefficients and equilibrium function); I, da, and B. 


Let Sr := L/TU, Ma := U /v^ and Re := UgL/p be the Strouhal, Mach and Reynolds 
numbers corresponding to hydrodynamics. Let also define C := c/U, L := Pda, Lg := d^/da, 
various dimensionless numbers corresponding to radiation. In all that follows, we assume our 
flow to be strongly under-relativistic so that C is large. 

Choosing B = I, we discover that the evolution of the dimensionless unknowns (still 
denoted in the same way) is governed by the following system of equations 

Sr dtg + div {gu) = 0, 

Srdtigu) + div{gu0u) + ;^(div (/rVM-k*Vn)+V(Adivn)) = 

' ^dtIo + ^dwh=Caa{B-Io), 

_ f dth + V/o = -C (da + Cgag) h, 

where g = g{t, x) G E+ and u = u{t, x) stand for the density and pressure, respectively, 
p = P{p) is the pressure, A = \{p) and p = p{p) are the viscosity coefficients. The given 
functions P, A and p are supposed sufficiently smooth, and we make the following strict 
ellipticity assumption 

V ■.= \ + 2p > d and p > d. 

In our recent work [8], we gave a mathematical justification of the low Mach number asymp¬ 
totics. In the present paper, we investigate another type of physically relevant asymptotic 
regimes, which are of diffusive type. They correspond to the case where C is large and all 
the other dimensionless numbers, but £ and Cg, are of order 1. To make it more concrete, 
take 

-1 


Ma = Sr = Re = 1, C = £ 


g = P'{g) = B'{g) = (^a{g) = crg{g) = i, 
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where e is a small positive number, bound to tend to 0. 

Because we shall focus on small perturbations of the reference density ^ = 1, it is conve¬ 
nient to introduce the new unknown h := B{q) — B{1). In this context, all the functions of 
Q may be written in terms of h. Setting jo := Iq — B{1) and ji := Ii, and using exponents 
to emphasize the dependency with respect to e, we eventually get the following system 

' dtb^ + + (1 + = 0 , 

dtiB + • VvF - (1 + k2{¥))AvF + (1 + k^{¥))V¥ = + ki{¥))% 

' edtil + H div jf = C{¥ - jg), 

edtjl + Vj§ = -£(1 Cs)% 

with A := /rA-|-(A-|-^)Vdiv and where /ci, k 2 , k^, k^ are smooth functions vanishing at 0. 


2. Formal asymptotics 

Let us first present some formal computations so as to exhibit the limit equations we can 
get from (|1.6p in different types of diffusive asymptotic regimes. We restrict to the case where 
the following necessary and sufficient linear stability condition (derived in [^) is fulfilled 

( 2 . 1 ) 

Note that (12.11) implies that liminf£e“^ > 0 for e going to 0. 

In all that follows, it is assumed that {¥, joAi) converges to (6,n, jo,ji) in some 
suitable space with enough regularity to pass to the limit in the nonlinear terms. 

2.1. Case A ~ e and Cg -(-oo. Denoting by V the orthogonal projector on divergence 
free vector fields, we get 

(2.2) Pj'!(i) = e-f(i+^=)*Pjf(0). 

Hence Vjf tends to () for e —0. 

2.1.1. Subcase CACg — )• 0. Setting Q := Id — P, we see that the equation for Jq entails that 
Qji = 0{e). Next, the equation for Qjl implies that Vjq goes to (), too, because e^Cg 0. 
Assuming that jo decays to 0 at infinity, this yields jo = 0. 

From the equation for jf, we also get 

(2.3) -C{l+Cg)]l = Vjl + 0{e). 

Hence e(l-|-iC<j)jf goes to 0 and (6, u) thus satisfies the barotropic Navier-Stokes equations. 
In other words, the radiative effect becomes negligible in the asymptotic £ ~ e and e^Cg —)• 0 
with Cg —)■ -t-oo 
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2.1.2. Subcase £ (0, +oo). This is the so-called nonequilibrium diffusion regime. 

The analysis of the previous paragraph shows that jf = 0{e) (hence ji = 0) and that (12.3h 
holds true. The new fact is that the equation for combined with (12.3p implies that 


(2.4) 




1 C 1 

n e C^Cs 


Ajo" = 0 { e ). 


Now, if we assume that 


C 


and 


C^Cs 


m 

9 ' 


e nu 

for some m G (0,-|-oo) and k>1 (see (|2.ip l. then (5, n) satisfies the following compressible 
Navier-Stokes equations coupled with a parabolic equation 


(2.5) 


dtb + u ■ V6 + (1 + A:i(6))div ft = 0, 

< dtu + u- Vft - (1 + k 2 {b))Au + (1 + k^{h))Vb + ^(1 + k 4 ,{h))V jq = 0, 
, dtjo + ^{3o-b-^/\jo) =0. 


2.1.3. Subcase C^Cg —)• +oo. We still have jf = 0{£), (j2.3p and thus (j2.4p holds true. Now, 
as C^Cs —>■ +00 and £ « e, the r.h.s. of (|2.4I) tends to 0. Therefore, if we assume as 
before that Cje Kjinv') then we find out that {b,u,jo) satishes the following degenerate 
nonequilibrium diffusion system 


( 2 . 6 ) 


dtb + u ■ Vb + {1 + ki{b))div u = 0, 

< dtu + u- Vu — (1 + k 2 {b))Au + (1 + k^{b))Vb + ^(1 + k 4 ,{b))V jq = 0, 
- %o + ;^(jo - 6) = 0. 


2.2. Case e <C £ 1. Recall that we have (12.3p while the equation for implies that 

(2.7) d\Y3l = nC{W - jff) + 0{£). 

Hence Qji = 0 (as £ —)• 0), and 

(2.8) AJq + n£^(l -|- Cs){b^ — Jq) = 0(e) 0(e£(l + £«))• 


Subcase CACg —)• 0. Then (|2.8p implies that Ajo = 0 and thus jo = 0 (if one assumes that 
jo —>• 0 at oo). Consequently, (12.3p implies that the radiative force in the velocity equation 
tends to 0 when e goes to 0. Therefore (6, ft) just satisfies the classical compressible Navier- 
Stokes equation. 


Subcase —>■ m G (0,-|-oo). We have ji = 0, and Relations ()2.3I) . ()2.8p imply that 

(6, ft, jo) fulfills the following Navier-Stokes-Poisson system 


(2.9) 


-|-ft • V6-|-(1-|-/ci(6))div ft = 0, 

< dtu^u- Vu- (1 -I- k2(b))Au -b (1 -b k^(h))Vb + ^(1 -b k4(b))Vjo = 0, 
, -i^^Ajo -b mn(jo - 6) = 0. 


Subcase C?Cg —)■ -boo. Then (j2.8p implies that jo = h. Combining with (12.31) . we thus find out 
that (b, u) fulfills the following compressible Navier-Stokes equation with modified pressure 
law 


( 2 . 10 ) 


dtb -b ft • V6 -b (1 -b A:i(6))div ft = 0, 

dtu + ft • Vfi - (1 -b k2{b))Au -b (l -b ^ -b k^(b) -b ^^4(^)) V6 = 0. 
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2.3. Case I ^ (0, +oo). 

Subcase —>■ m E [0, +oo). Passing to the limit in (12.8p gives 

(2.11) - + n(£^+ m)(jo - &) = 0. 

So we get System (I2.9jl for {b,jo,u) with the last equation replaced by (12.lip . 

Subcase Cg —)• +oo. Exactly as in the case /I —>■ 0, we get jo = b, ji = 0, and (6, u) satisfies 

(ICT) . 

2.4. Case C +oo. Relation (12.3p implies that ji = 0, and thus, according to (12.7p . we 
have jo = b. Therefore (12.3p implies that 

£(l+£,)jf ^V6, 

and (6, n) thus satisfies (|2.10l) . 

To make a long story short, the above formal computations pointed out five types of 
asymptotic regimes. They are governed by 

(1) The ordinary compressible Navier-Stokes equations with null radiation (if £ —)• 0 and 

Cg —>-0); 

(2) The compressible Navier-Stokes equation with an extra pressure term see (|2.10l) (equi¬ 
librium diffusion regime corresponding to e <C £ and C^Cg —>■ -|-oo, or £ —)• -|-oo); 

(3) The Navier-Stokes-Poisson equations (12.9p (or (j2.1ip l (case e <C £ < 1 and v‘^C?C.g 
m E (0, -boo)); 

(4) The compressible Navier-Stokes equations coupled with a parabolic equation (j2.5p 
(nonequilibrium diffusion regime £ « e and v^CgC^ —>■ m E (0, -boo)); 

(5) The compressible Navier-Stokes equations coupled with a damped equation (j2.6p 
(degenerate nonequilibrium diffusion regime £ ps e and CgC'^ —)• -boo). 

The rest of the paper is devoted to justifying rigorously the last four asymptotics globally 
in time in the framework of small solutions with critical regularity. In the next section, we 
introduce a few notations that will be needed to define our functional framework, and give 
an overview of the strategy. Section 0] is devoted to a fine analysis of the linearized equations 
(11.61) about (0,0,0,0), which turns out to be essentially the key to proving global results and 
justifying the diffusive asymptotics we have in mind. The next three sections are devoted to 
the rigorous justification of the nonequilibrium diffusion regime £ « e and CgC^ > 1, the 
equilibrium diffusion regime £ —-boo and of the Poisson type diffusion regime (e <C £ < 1 
and Hg —)• m E (0, -boo)). In all of those sections, we establish a global-in-time existence 
result for the expected limit system, and for (11.61) supplemented with uniform estimates (for 
coefficients £ and Cg satisfying the assumptions of the studied regime), and eventually 
show the convergence of the solutions of (jl.6j) to those of the expected limit system. Some 
estimates, of independent interest, for the solutions to a class of linear ODEs corresponding 
to the linearized equations of (II.6j) in the Fourier space are postponed in the appendix. 

3. Functional framework and overview of the method 

The functional framework we shall work in is modeled on the linearized equations cor¬ 
responding to (11.61) . and is thus the same as in our first paper [7] devoted to the global 
well-posedness issue in critical regularity spaces for small perturbations of a stable constant 
state. The key to proving asymptotic results however, is to prescribe norms depending on 
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the parameters e, £ and Cg, so as to get optimal uniform estimates, enabling our justifying 
rigorously the different diffusive asymptotics exhibited above. 

Let us first very briefly recall the definition of homogeneous Besov spaces ^ (the reader 
is referred to [T], Chap. 2 for more details). For simplicity, we focus on the M"" case (adapting 
the construction to the torus being quite straightforward). Fix some smooth radial bump 
function x : ^ [ 0 , 1 ] with x = 1 on i?( 0 , 1 / 2 ) and X = 0 outside B( 0 , 1 ), nonincreasing 

with respect to the radial variable. Let ip(X) := x(?/2) — x(0- The elementary spectral 
cut-off operator entering in the Littlewood-Paley decomposition is defined by 

Aju := (p{2-W)u = T-^{ip{2-W)Tu), j E Z 

where we denote by J- the standard Fourier transform in 

For any s E M, the homogeneous Besov space i is the set of tempered distributions u 
so that 

||tt||gs := ^ 2-^*||Ajn|A2 < oo, 
jez 

and 

(3.1) lim x{^D)u = 0 in L°°. 

A—^“)~CXD 

As pointed out in [7] , scaling considerations that neglect low order terms of System (jl. 6 p 

, n_^ . n 

suggest that critical regularity is B^i for uq, jq^o and ji^o, and B 21 for 6 o- However, to 
handle lower order terms, one has to make additional assumptions for the low frequencies. To 
this end, it is convenient to introduce the following notation (where r/ stands for a positive 
parameter) 

||u|||/ := 2 '=^||A,u|U 2 and ||u||^’/ := 2 '=^||A,u|U 2 , 

^2,1 f —' ^2,1 ,—' 

2>‘<2r) 2^>ril2 

and also 

:= ^ AfcU and := ^ A^u. 

2*‘<rj 2^>r} 

Note that ||u^’^||rs < C\\u\f^ and ||u^’^||rs < . Because the Littlewood-Paley 

"-“24 "-° 2,1 

decomposition is not quite orthogonal, it is important to allow for a small overlap in the 
above definition of norms. 

In some places, we will have to specify also the behavior for the middle frequencies, by 
considering for given 0 < rj < r]', 

r]<2^<r]' 

Broadly speaking, our strategy to justify the different types of diffusive limits is as follows: 

• Step 1: We prove ‘uniform estimates’ for the global solutions to (|1.6I) . uniform mean¬ 
ing that we want a bound independent of e, but the norm itself may depend ‘in a 
nice way’ of the parameters e, L and Ca- 

• Step 2 : We show that the limit system is globally well-posed in the small data case. 

• Step 3: We take advantage of estimates of Step 1 to exhibit weak compactness prop¬ 
erties. Combining with the uniqueness result of Step 2 , this allows to conclude to the 
convergence of the whole family of solutions of (ll. 6 jl to those of the limit system. 
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The most technical part is step 1, as it requires a fine analysis of the linearized equations of 
(jl.6p about 0 that keeps track of the coefficients £, and e. Schematically, in the Fourier 
space, one has to resort to different types of estimates for low, medium and high frequen¬ 
cies. The low frequency analysis is carried out by considering approximate eigenmodes of 
the system, that are constructed by a perturbative method from the (explicit) eigenmodes 
corresponding to null frequency. A part of the difficulty is that the ‘fluid modes’ are of para¬ 
bolic type, hence the corresponding eigenvalues tend quadratically to 0 when the frequency 
size tends to 0 while the radiative modes are expected to be exponentially damped. The 
high frequency analysis is inspired by the corresponding one for the barotropic Navier-Stokes 
equations, after noticing that coupling between radiative and fluid unknowns occurs only 
through 0 order terms, and thus tend to be negligible for very high frequencies. Last but not 
least, medium frequency regime has to be looked at with the greatest care, as the low and 
high frequency regimes need not overlap. There is no general strategy for handling them, 
apart from guessing approximate eigenmodes of the system. 


4. Uniform estimates for the linearized equations 

In order to reduce the study to the case where the total viscosity u := \ +2^ is 1, and to 
get a symmetric first order system for the radiative unknowns, let us set 

(4.1) {b,u,jo,ji){t,x) := (6^u^ v^jo, jf)(vLi/x). 

Then , Jq, jf) satisfies (|1.6p if and only if {b,u, jo, ji) satisfies 

dtb + u ■ Vb + {1 + ki{b))div u = 0, 

dtu + u-Vu - {1 + k 2 {b))Au + (1 + /c3(6))V6 = %^(1 + ^4(6))^!, 

£dtjo += C{b - y/n jo), 

edji + = -CMji, 

with 


(4.3) 


M.-.= 1 + Cs, C:=i>C and A:=i' ^A. 


The corresponding linearized system reads 

dtb + divu = /, 


(4.4) 


dtu -Au + Vb- ^ji = g, 

£dtjo + ^ div ji -k C{jo - ^/nb) = 0, 
edji -k ^ Vjo + CMji = 0. 


On one hand, the coupling between the incompressible part of u and ji that is Vu and Vji 
where V stands for the projector on divergence-free vector-fields is obvious as 


(4.5) dtPu— —AVu = - Vji, 

V n 

and 

'Pji{t) = e e Vji{0), 
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hence in any functional space X we have 

(4.6) CM\\Vj4LHx)<e\\Vjim\x. 

On the other hand, the coupling between 6 , d := A“^divn, jo and ji := A“^divji (where 
A^ := (—is quite complicated: in Fourier variables, we have 


(4.7) 


A 

dt 


(^) 


/ 

0 

P 

0 

0 

d 



-p 

P^ 

0 

CM 

n 

+ 


y/nC 

0 

c 

p 

Jo 



£ 

£ 

6^/n 

V ji y 


V 

0 

0 

p 

CM 


ey/n 

£ 



( ^ ^ 




d 


0 


Jo 


0 


\3i J 


VO/ 


The analysis that has been performed in [7] pointed out the following necessary and sufficient 
stability condition 

(4.8) C> -{1 + M-^). 

n 

So we shall make this assumption in all that follows. Of course one also has to keep in mind 
that A4 > 1, a consequence of Ai := l + Cg- For notational simplicity, we shall simply denote 
A by £ in the following computations. 


4.1. Estimates for small p. In order to prove estimates in the case 0 < p < Ci (with 
Cl > \/l + ), we shall use that (I4.7jl enters in the class of ODEs that has been considered 

in the Appendix. Indeed, it corresponds to (IA.3p with 


(4.9) 


CM 


? = 


n 


rj = 


y/nC 


/? = 


C 


1 


a = 


evra 


7 = 


CM 


4.1.1. The case C > 1 and CeM > 1. We shall follow the first approach proposed in 
Appendix O It corresponds to the following matrices Aq, 


An — 


( 0 
0 
0 
V 0 


0 

0 
c 
£ 

0 0 


Bi = - 


( 0 
0 
0 

Vi 


0 

0 

0 

CM 


\ 


0 

0 


n 


0 

1 

1372 

0 

0 


/ 


= 


- \ 
n \ 

0 
0 
0 / 


0 

-1 - 

0 

0 


and A 2 = 


Ai, 

A 2 and Bi 


1 

0 

0 

\ 

0 

0 

0 


0 

0 

l+£= 

y/ne 


0 

1 

'Jn£ 

0 

/ 

/ 0 

0 0 

0 V 


0 

1 0 

_ £_ 

n 


0 

0 0 

0 


VO 

0 0 

0 y 



Therefore we set 


fA in'! 


p — 


0 

0 


0 


nCM 

0 
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which corresponds to the change of unknowns 

(4.11) 


(^ ^ 


/ 

1 

0 

0 

£2 

uCmP 

\ 


( ^ ^ 

b 



-^p 

1 


£_ 

n 



d 

jo 



-\/n 

-^0 
£ P 

1 

£2 

V^cP 



k 

V / 


V 

~'CmP 

0 

0 

1 

/ 


\ A / 


According to ()A. 8 h . working with (a, d, jo, ji) or (b, is equivalent whenever 

(4.12) p < £min(e“^, A4). 

Let us first compute the matrices PBi, [P, Ai] and A 3 appearing in (IA.2P 

/ -M-^ 0 0 0 \ 


nC 


V 


0 

0 

0 


-10 0 
0 1 0 

0 0 M-^ ) 




£ 


0 0 

0 0 

EyJnM. 


= 


n 


0 


£ 

CM 


-l-Af-i 


V 


(l+n-l)£2 

0 

0 


e 

TP 

0 


£( 1 +X- 1 ) 

„ 3/2 

0 

0 

0 

0 


Cs/n £2^3/2 

0 0 
n _^ 


l+£2(l+AI-l(n+l)) 

-- 

0 

0 

\ 

nTPJTP 

0 

£( 1 +£ 2 ) 

0 


/ 


Because £ > 1, we thus have IA 3 I < Hence, up to a 0{ep^ jC) term, the equations for 
(b,c)) read 

d ' 


dt 


+ P 


0 1 
-l-n"^ 0 




= P 


nCM 
0 

£(1+A4-1) 

n^/^C 

0 


0 

1 _^ 

nC 


£ _ 1+£^(1+A4 ^(n+1)) 
n rPT 


jo 


In order to estimate (b, b), we just follow the method of Appendix[Bl which requires Condition 
(14.8p and 

Vl + n~^ 


(4.13) 


P< 


^ nC 


Keeping (I4.12h in mind and noticing that 


P = 1-4(1+M-') 
nL 

is of order 1 for small e, we thus conclude that if 
(4.14) 


p < \/\ + n ^ and / 9 <£min(e ^,A4), 
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then 


(4.15) \ibM)\+p" / |(b,h)|dr<|(b,0)(0)| 


+ P‘ 


I(i“ 


+ |ji| ) dr + 


££_ 

C 


,3 ft _ 

r |(b,c)Jo,ji)| 


dr. 


Next, we see that the equations for (jo,ji) read (omitting the 0{£p^jC) term) 

^2 


(4.16) 


d 


io 


dt I ji 


+ 


C I ep" 
e nC 

0 


0 

CM I 
e nCM 


1 ° 


+ 


P 


n e 


0 1 + e^ 


-1 


= P 


l+£^(l+At(l+n)) 

e^/nCM 

0 


0 

0 

1+M' 


jlo 

ii 

b 

b 


Therefore, computing 

(4.17) 


d 

dt 


liol^ + (1 + ^^)|ji| 


so as to eliminate the term in p, we end up with 

(4.18) |0o,ji)(t)| + - / |io,ji|dr < |Oo,ji)(0)| 

£ Jo 


+ 




Now, adding up (I4.15P and (|4.18l) . we easily conclude that if e is small enough and 

(4.19) £min(l,eA4) > 1, 

then we have 

£ 7* 

e Jo 


(4.20) |(b,'D,jo,ji)(t)|+/ 0 ^ / |(b, 0 )|dr + 

Jo 

whenever 0 < p < y/l + n~^. 


IOo,ii)|dT < |(b,t),jo,ii)( 0 )|. 


Now, resuming to the ji equation in (j4.16p . and evaluating the hrst order term according 
to (I4.20p . we deduce that, in addition 


(4.21) 


CM 


\h\dT < |(b,c),jo,ii)(0)|. 


4.1.2. The case e <C £ < 1 with eC^C}^ 1 and C'^Cg 1. If £ <C 1 then plugging (|4.15p 
in (j4.18p does not allow to get (I4.20p any longer. In order to overcome this, we shall follow 
the second approach proposed in Appendix O with coefficients defined as in (14.9p : we set 


P = 


( 0 
0 


0 


CM 


0 

0 

epn 


0 

0 

1 

y/n(l—M)C 


nCM 

0 

l+£^ 


\ 


^/n(l—M)C 

0 
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and we thus have, remembering that Ai — 1 = Cs 

( b\ / 1 0 


(4.22) 


V = 


d 
jo 

vTi y 


nCf 


1 




\ CMCs 

The determinant of the above matrix is 


, \/ne 1 

-\/n 1 

0 


y/nCsC 


uCmP 

£ 

n 

AnCCs ^ 

1 


\ / n 

d 

h 

/ V y 


1 + 




n £2 


P 


1 + 




1 + e^ 

nC^C^J 


Hence working with {b,d,jo,ji) or (b,t),jo,ji) is equivalent whenever 


(4.23) 


C 

p < — and p < y/n CCg 

£ 


Then following the computations of Appendix!^ second approach, and setting A 3 ;= (PAq — 
Ai)P^ + A 2 P leads to 


(4.24) -V + p 


+ 


= P 


( nCM P^ 

0 ( 1 - 

0 

0 

0 

0 

l+£^ _ (n+l)£ 

e^paCaC y/nC 

0 


1 

0 

0 

\ 

0 

0 

0 


0 

0 

0 


0 

0 

0 

y 

0 




e 

nC 

) p ^ 




V 


0 

0 




l+£^ 


0 

0 

0 

1 


CCaM 


{l+M-pe 
rpl'^C. 

Q A _ !+£• 

0 
0 


nCM neCC: 

0 \ 

_^ _ £^(l+n) 

n n^C TTCAT 

0 


■)p^) 


V 


0 


V 




+ p\l + pP)A^{I + pP)-^V. 


Let us bound A 3 in order to determinate for which values of p the last term in (j4.24p is 
indeed negligible. Just writing that IA 3 I < |Pp(|P| |Ao| + |Ai|) + IA 2 I |P| using the explicit 
values of Aq, Ai and A 2 and 


(4.25) 


^1 ^ 4 max( e, ^ 


C 


\ T 

V •‘-'i 


does not provide an accurate enough bound for A 3 . Hence one has to go to further compu¬ 
tations. Now, we get 


/ 


PAnP^ = 


0 

l+£^ 

n^C^MCa 

0 


£( 1 +P) 

\ nC'^MC'i 


£" 


0 

+£■ 

PJEcPEl 


( 1 +£^)A< 


0 

1 { l+£^ _ p2\ 

0 




£ 

M 


') \ 


0 


0 


Pl'^C.aC? 


£ - 


(l+e^)M ( £2 _ 1+£2 \ 
rPI'^eCsC? VAH ~TP> 

tt) “ 
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AoP = 


/ 0 0 0 0 \ 

£ n £ n 

nCM ^ 


0 0 0 0 

\ 0 0 0 0 / 


f 0 

(l+n-l)£2 

0 

V 0 


£2 _£( 1 +£!)\ 

nC'^ nC-^C. 

(l+n“l)e 2 

n^/^C^CsM 

0 0 


0 


and AiP^ = 

0 

0 y 

Hence, given that e <C £ < 1 and that Cg M in the regime that we are considering, one 
may conclude that 

el 1 


( 4 . 26 ) 


IHsI < max 


CALg eC?C?g, 

Note that we still have u ^ 1 for e ^ 0. Hence applying the method of the appendix to 
handle (b, t)), we find out that if (|4.8I) . (|4.13l) and (|4.23l) are fulfilled then 


i(b,0)(t)i+p2 ri(b,b)idT<i(b,0)(o)i+p2 

4o 


As regards the radiative modes, we have 


^liol + ^Ihl ]dT 


lio(i)| + 


c 


+ (e + 


i±£!^Z 

eCg / nC, 


\]o\dT < |jo( 0 )| pCp^ 

+Cp^ max 


i:u 




dr. 


eCCs ^ £ 
1 


dr 


C C?Lg eC? a 


yyVirfr, 


liiWI + 


PM 


+ 


1 +£2 


nCM neCC, 


P 


\h\dT < |ji( 0 )| +C 


P 


CCsM 
£ 1 


| 0 | dr 


+Cp3m»c(i Vldr. 


From the above three inequalities, we get for any A G (0,1] 


|(b, b)(t)| + A|jo(t)| + |ii(i)| + P‘ 


r 1^5)1 

Jo 


dr + A— f liol dr + 
£ Jo 


CCg 


e Jo 


\]i\dT 


|(b, t))( 0 )| + A|jo( 0 )| + |ji( 0 )| + ^ Jq (eZc” z) 


+ 


r 

mJo 


11 )| dr + (? max 


1 


£’ £2£ ’ e£2£2 




dr. 


Now, we notice that taking A = cq min(1, £ 12 / 15 ) for a sufficiently small constant cq allows 
to absorb all the terms of the r.h.s. (but the data) by the l.h.s. provided we have 0 < p < 
\/l + n~^ 

(4.27) £ < /: < 1, ££ 2£2 ^ I ^ £ 

We thus conclude that for all 0 < p < Vl + n~^, we have 

(4.28) |(b,b,min(l,£££ 5 )jo,ji)(t)| + p^ / |(b, 0 )| dr + — min(l, 

do ^ 


eCCg 


Do I dr 


+ 


CM 


Dll dr < C|(b,b,min(l,££/l 5 )jo,ji)( 0 )|. 
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Let US point out that in the case where C^Cg ~ 1 (even if ~ e in fact) then the same 
computation will lead to (|4.28p . but only for 0 < p < c, with c a small enough constant. 

4.1.3. The case e <C and C^Cg ~ 1. As the value y/l + n~^ will not play any 

particular role, we fix some Ci > c in the following computations. We want to get (j4.28p for 
p G [c, Cl]. To this end, we introduce Co aiid Ci such that 


Co jo - 


n - ^ ^ p ^ 

i + sfer 


Then we discover that (6, d, CoiCi) fulfills 

d$+ pd = 0, 

dtd + p^d - pi I + 


1 


n 


dtCo H-Co =-^Ci - Vn -^^—2— d, 

e Ey/n 1 - P. 


^ Co + ^Ci, 


3/2 


n 


P 


dtCi + 


fv- 


M- 


nC?M 


TcTm 


Ci = 


1 + 


yJneM. V nCT'M. 


Co- 


Let p be in [c, Cij. For the first two equations, performing the standard barotropic estimates 
(which rely on the use of Up defined in (14.431) ) leads to 


(4.29) |(6,d)(t)|+/ |(6,d)|dr<|(6,d)(0)|+ / jColdT + CC, / jCildr. 

Jo Jo Jo 

For Co) it is obvious that 

(4.30) \Co{t)\ + - f IColdr < |Co(0)| +-^ / \Ci\dT + py/nf \d\dT, 

e Jo Jo Jo 

and, as our conditions on U and Tg guarantee that p < y/rij2CJA for small enough e, we 
also have 


(4.31) 


ICi(t)l + 


CM 


Ticii 

Jo 


dr < |Ci(0)| + 


P 


y/neM 




jdj dr. 


Putting together those three inequalities, we readily get for all A,B > 0, observing that 
PS C‘^M ^ 1 


\{b,d){t)\ + A^|Co(t)| + Be\Ci{t)\ + \(b,d)\dT 


rt 


rt 


+ A / \Co\dT + BCM / ICijdT 


< |(6,d)(0)|+A-|Co(0)|+i?e|Ci(0)|+CM / jCil^r 


+ 


^ + \Ci\dr+ A^J^ |dldr + ^^ ICojdr. 


It is now clear that if one takes first A large enough (independently of e) and B much larger, 
all the integrals of the r.h.s. may be absorbed by the left-hand side, as M~^ <C p, and we 
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thus get for all p E [c, Ci] 

(4.32) |(6,d)(t)| + ^|Co(i)|+e|Ci(i)l + ^ \(b,d)\dT 

+ [ \C,o\dT + CM f ICil^ |(&,c?)(0)| + —|Co(0)| + £|Ci(0)|- 
Jo Jo ^ 

Plugging this new inequality in (j4.3ip . we easily deduce that 

ICi(t)l + ^J‘ Ifil 

then inserting this information and (14.321) in (|4.30p . we discover that 

I Co (0)1 + — ^ I Co I dp ^ I Co (0)1 + £|Ci(0)| + 

Therefore, putting (I4.32p and the above two inequalities together, using that |(&, d, Co) Ci)l ~ 
1(6, d, jo, ji)| and assuming in addition that C < we conclude that 

_/■*_ C cc /■* -- 

(4.33) |(6,d,-jo,ji)(t)| + / |(6,d)|dr + - / |Co|drH- -\ |Ci|dr 

^ Jo ^ Jo ^ Jo 

^ \{b, d, ^jo, ji)(0)| for all c < p < Ci. 


Note that due to the expression of Ci, one may replace ("i with ji of ji in the integral. 
Still in the case C < we claim that we have the following inequality 


(4.34) 


\pjo{t)\ + 


C 


IpCoI dr < |(6,d,pjo, ^jo, ji)(0)| for all 0 < p < Ci, 


which turns out to be crucial in the justification of the asymptotics toward (j2.9p . 

Indeed, inequality (|4.30p does not require any assumption on p and thus implies that 

C 

e 


f* pt 2 pt pt 

\pCo{t)\ + - \pCo\dT <\pCo{ 0 )\ +\Cl\dT + p^^/n \d\ 

£ Jo Jo Jo 


dr. 


For 0 < p < c (resp. c < p < Ci), the last term may be bounded according to (I4.28h (resp. 
(j4.33p l. Regarding the term with (^i, we notice that 

hence (I4.28P and the fact that C^Cg ~ 1 guarantee that for 0 < p < c 


eWn 


ICildr < C 


CM 


liil dr + 


£2 


P^\jo + '/nb\ dr 


< 


£(l + ^)|(6,d,e££,jo,ji)(0)|. 


In the case C < it is obvious that Inequality (|4.33l) implies that the above inequality is 
also true in the range c < p < Ci, which completes the proof of (I4.34h . 
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4.1.4. The case £ « e and CsS^ > 1. We saw that if (I4.8p is fulfilled then ()4.28ll holds true 
on some small interval [0, c]. So we have to fill in the gap between c and Vl + n~^. As the 
value Vl + n~^ does not play any particular role, we fix some C\> c and look for estimates 
if /? G [c, Cl]. For simplicity, take L = Kejn (with k > 1 owing to (|4.8p l. 

Setting Cl := ji - ^CM ^ ~ as before, we observe that 


(4.35) 


dtb + pd = 0, 

dtd + p‘^d- pb = ^jo + ^Ci, 
dtjo + (f + + ;;|^Ci = -^b, 

, dtCl + ) Cl = K,'^e-‘M)dO ~ 

Let us focus on the subsystem corresponding to the first three equations, namely 

dtb + pd = 0, 

did + p^d- pb - ^jo = /, 

^ dtjo + (^ + = g- 

If we have the stronger condition e^A4 —>■ oo then we rewrite System (|4.36p as follows 


(4.36) 


(4.37) 


jL 

dt 



, / 


+ 

\jo) 

\ 


—^ 0 

Vn. 



0 
/ 

p2 ^ ^ 

'TF^do +9J 


The eigenvalues of the matrix Mp in the left-hand side are the roots of the polynomial 
-I- ai(p)A^ - 02 (p) A -I- 03 (p) with 

al(/5) = -+P^ a2(p) = (l +-)p^ 03 (p) = fl +-)—• 

n \ nJ \ nJ n 

According to Routh-Hurwitz criterion, those roots have positive real part if and only if 


ai(p) > 0, 


ai(p) 1 
03(p) 02(p) 


> 0 and 


oi(p) 1 0 

03 (p) 02 (p) oi(p) 

0 0 03 (p) 


> 0 . 


As oi(p) and 03 (p) are positive, it suffices to check the second condition, that is 

oi(p)o2(p) - 03 (p) = (l + -)p‘^ + - 1) > 0, 

V n/ 

and this is indeed the case for all p > 0, as k > 1 . 


In particular, all the eigenvalues of the matrix Mp have positive real part if we assume p 
to belong to the compact set [c, Cij. Therefore (see [7]) there exist two positive constants C 2 
and C 2 depending only on c and Ci, so that the matrix Mp satisfies 

\^-tMp\ < for all t > 0 and p G [c,Ci]. 

By taking advantage of Duhamel’s formula, we thus deduce that 

|(6,d,j/o)(t)||(a,d,jo)|dT < |(6,(i,jo)(0 )|^ j^ IjoIc^'T- 
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Of course, owing to the assumption oo, the last term of the r.h.s. may be absorbed 

by the l.h.s., for e going to 0. So we get 

(4.38) l(b,d,jo)(t)l + f l(b,d,jo)ldT<l(b,d,jo)(0}l+f l(f,g)ldT. 

Jo Jo 

In the case where does not go to oo then we have to proceed slightly differently. If we 

assume (for simplicity) that e^A4 tends to some m > 0, then the matrix Mp in (j4.37p has 
to be changed in 

'0/9 0 


Np = 


-P 

^ Q — 

y/n n ' Km, 


The above analysis based on Routh-Hurwitz theorem still holds as the additional term has 
‘the good sign’, and one may conclude, as before, that (j4.38l) is satisfied for all p G [c, Ci]. 

In every case > 1, resuming to (j4.35p . Inequality (|4.38p allows us to get for all 

P G [c,Ci] 

(4.39) I(6, d, jo)(i)I + / \(J>,d, jo) \ dr < \ (b, d, jo)(0)| + eA4 / |CiI 

Jo Jo 

Next, from the equation of Ci, we readily get 

\Ci{t)\+J\4 f \Ci\ dr < \Ci{0)\ + f \(b,jo)\dT. 

Jo Jo 

Hence, adding up to Inequality (14.391) . we conclude that (I4.28P is also true for all p G [c, Ci]. 
This completes the proof of estimates in the low frequency regime p G [0,(71]. 

4.2. Estimates for middle frequencies. 

The case liminf A4 = +oo. As in the previous paragraph, introduce ("i := ji ~ ' 

The system fulfilled by (6, d, jo,Ci) reads 
dtb + pd = 0, 
dtd + p^d- pb = ^jo + 
dtjo + i (^ + T^)jo + ^Ci = y/njb, 

^ dtCl + ^ ~ 

The subsystem corresponding to the first three equations is 

dtb + pd = 0, 

did + p‘^d- pb - ^ jo = /, 

. 5iio + 7(l + ;^)jo-= 


(4.40) 


(4.41) 


with 

(4.42) 


/ =-Cl and g = - ^Ci- 


n 


n e 


Assume that (/,p) = (0,0) for a while and set 


(4.43) Ul := 2|(6,d)|2 - 2pRe(6d) + \pbf 
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On one hand, we have 


(4.44) 

^ ^ 2dt 

and on the other hand, 

1 d 


+ p‘^\W + p‘^\d\^ = {{‘^d - pb) jo), 


C 




,^2^;boP = 77fRe(6jo). 


Therefore 




= 2^^Re (d Jo)- 


Now, by using the fact that 




2 , ^p r i 2 


and by taking A = 3/4, we conclude that for > 3 ( 4 ^^yyj we have 

+ £2l(^>^oo)P < 0, 

whence, because ^{An^-i) — ^ n’ some universal positive constants cq and C 

\{pb,d){t)\ + ^J~^\pjo{t)\ < Ce~^°^(^\{pb,d){0)\ + ^J^\pjoiO)\^ for p > \/l + n-^- 

Resuming to the equation fulfilled by jo in (I4.4ip . the above inequality implies (still assuming 
that f = g = 0) that 


^11 + 


I’ ' ) 

nC?M) 


\pjo\ dr < |pjo(0)l + 


'0 


y/nC /■* 

e Jo 


\pb\ dr 


< W- 


\pjo{0)\ + ^l(d,p6)(0)l. 


then plugging this inequality in the equation for d, we get in addition 

rt 


P^ [ ^ |rf( 0 )| + |/3%)| 

Jo 


+ \/^l/5jo(0)|. 


Repeating the above computations in the case of general source terms / and g, we conclude 
that the solution (5, d, jo) to (I4.4ip satisfies for all p > Vl + n~^, assuming only that C> e 


\{pb,d){t)\ + ^\pjo{t)\ + \{pb,p^d)\dT + j 


1 + 


nC'^M 


\pjo\ dr 


< C|^|(p6,d)(0)| + y ^|pjo(0)| + (^l/l + y ) dr 

Resuming to the value of / and g in ^4.421) . we thus get for p > Vl + n~^ and £ > 


(4.45) |(p6,d)(t)| + y^|pjo(i)| + \{pb, p‘^d)\dr 




nC^M J Jo 


< C ( \{pb, d)(0)| + W w|/3jo(0)| + ( £A4 + 


1/5 jo I dr 


\fCe 


\Ci\ dr 
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Next, it is clear that the equation for (^i implies that whenever p < CM 


(4.46) |Ci(t)| + 


CM 


Ticii 

Jo 


dr < |Ci(0)| + - 


1 + 


y/neM \ nC^M ) Jq 


[ \pjo\ 
Jo 


dT + 


eM 


\ph\ dr. 


Hence, we get if M is large enough and p^ <C C^l‘^ 


(4.47) 


ICi(i)l + 


CM 


|Ci|^^r<C |Ci(0)| + 


eM 


pb,d){0)\ + \pjo{0 


Then plugging that inequality in (14.451) implies that 

(4.48) \{pb,d){t)\ + ^\pjo{t)\+p‘^J^\d\dT + J^ \pb\dT 




/Viol 

40 


dr 


nC^Mj Jo 

< c(|{Arf)(0)| + y|lpio{0)| + (^ + ^^)lfi(0)l 


whenever l + l/n</9^<C 

Here is another method that gives decay estimates in the range C^fM <C p <C CM if M 
is large enough. From the first two equations of (|4.40p . we have 

(4.49) |(p6,d)(t)|+p^ / |dldT+ / \pb\dT <\{pb,d){0)\ + -^ [ \jo\dT + ^J^ [ \Ci\dT. 

Jo Jo Jo n Jo 

The equations for jo and Ci give, if /? < CM 


(4.50) |jo(^)| H (4^ + 


P 


nCMJ Jo 


\jo\ dr < |jo(0)| + 


P 


(4.51) |Ci(t)| + 


CM 


iCil dr < |Ci(0)| + 


P 


C + - 


ne Jo 
P^ 


Ticii 

4o 


dr + 


Cy/n 


I Jo I dT + 


eM 


\b\ dr, 


|6| dr. 


y/n CMe \ nCM ) Jq 
Plugging (I4.50p in (j4.5ip , we discover if p <C CM that 

(4.52) |C,(«)| + ^l‘&\dT < |fi(0)| + ^\j,{0)\ + ^l\\b\dT. 

Inserting that inequality in (I4.49h , we conclude that the last term of (j4.49p may be absorbed 
by the l.h.s. if M is large enough. Now, Inequality (I4.50p guarantees that 

f p^CM WVij- ( ^ '' 

^\n^Cm + p^n) Jo \nCm + p\ 

Again, resuming to (14.491) . we see that the second term in the r.h.s. may be absorbed by the 
l.h.s. This is also the case of the last one if p^ ^ C?‘M. If all those conditions are fulfilled 
then we end up if C^J~M <C p < CM with 


n 


3/2 


\\loV 

Jo 


t p\b\ 

Jo 


dr. 


(4.53) \{pb,d){t)\ + p^ j \d\dT + pl \b\dT + p 


/ |6| dr + p / I 
Jo Jo 


dr 


^ _ eCM - - 

+ CM / |Ci|dr<|(p6,d)(0)| + —-|jo(0)|+e|Ci(0)|. 
Jo P 
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Resuming to (14.501) and (I4.52p . we thus easily deduce first that 


and next that 


ICi(^)l + 


CM 


Ticii 

4o 


dr < |Ci( 0 )| + bo( 0 )| + ^|(p&,d)( 0 )|, 


bo(OI + “ (-^ + 


p 


CM 


lioUr < |j„(0)| +^|Ci(0)| + i|(p6,<i)(0)|. 


Our estimates and the definition of allow us to change to ji in (|4.53l) . So finally, in 
the case eM ^ 1 we get for some large enough Ci and small enough c independent of e 

(4.54) \{pb,d,{l + ^^)jo,ji){t)\+ [ IpblclT + p"^ f \d\dT 

Jo Jo 


+ 


r(l + ^) 


/ \pjo\ 
Jo 


dr + 


CM 


eCM ^ ' p 

whenever CiC'/M < p < cCM. 

The case e^M <1/2 and £ S> e. Let Co := Jo ~ ^/nb. We start from 

dtb + pd = 0, 
dtd + p^d- pb = ^Ji, 
dtCo + fCo + - y/npd = 0 , 

^ 9ji + ^Ii-^Co-f 6 = 0 . 

In order to show the exponential decay, we set 


ICil dr < \{pbX (l + ^)jo,Ji)(0)|, 


(4.55) 


Ul := 2|6|2 + 2|d|2 + \pb\^ - 2pRe (bd) and := |CoP + |jip. 


We easily get for all K >0, 


(4.56) 


2 dt 


C 


ui + iLJ-/ + p^\{b, d)\^ + K- iCor + -M|J 1 


= 2 (^^Re (ji d) + p(— - ^^^^Re {jib) + K^/nRe {pdCo)- 
n \ £ n J 

It is thus natural to take K = ^CM to cancel out the second term of the r.h.s. For the first 
and the last terms, we write that 

Ky/nRe (pdCo) < \p^\d\‘^ + nK‘^\jo - ^/nb\‘^, 

^Re(jid) = ^Re(U) < y^\d\^ + 

Note that the last terms above may be absorbed by the l.h.s. of (|4.56l) if, say 

r 2K 

nK < — and -< CMf?. 

- ‘le e - ^ 

Given the value of K, the first condition is equivalent to e'^M < 1/2, whereas the second 
one means that p^n > 2. Under this latter condition, we thus end up with 


(4.57) 




eCM 


n 


1 


J^^)+-p^\{b,d)\^ + 


C^M 


n 


ICol' + xA^Ijin <0, 


























20 


RAPHAEL DANCHIN* , BERNARD DUCOMET** 


which implies, according to (|4.8p . the following exponential decay estimate for some small 
enough k > 0 


(4.58) Ul{t) + eCMJlit) < {uHq) + 


eCM 


n 


Jo 


if p > -s/Jjn. 


To exhibit the parabolic decay for d, we introduce Ci := Ji — [y/nCM.) ^pjo, and get 

dih + pd = 0 , 

dtd + p^d - p{l + ^)b = + ;^Co, 

dtC,o + = y/npd - 

, dtCl + ^ y/^CMe “*■ TSV()^0 + neC'^M^^- 

We thus have for yfJJn < p < \pnj2 CM. 


r ip^i 

Jo 


dr + p"^ \d\dT <\{pb,d){0)\ + 


^ f \Ci\dr+ ^ f\Co\dr, 


-ic + 


nCM J Jq 
CM 


ICol dr < |Co(0)| + 


y/ne 


iCil dr + y/np / |d| dr + 


y/neCM Jo 


\pb\ dr, 


2e 


ICil dr < |Ci(0)| + 


y/neCM 


C + 


P 


nCM J Jo 


ICol dr + 


P 


neC^M? 


\pb\ dr. 


Combining the inequalities for and (Ci, we easily get ii p < cCM with c small enough 


( ICol fir <e|Co(0)|+ -£^|Ci(0)1+pe^ Idjdr 


CM 


+ 


[ \pb\c 
Jo 


dr. 


Now, the exponential decay pointed out in (|4.58p allows to bound the last terms above, and 
we get 


(4.59) 

(4.60) 


/ iCildr <^/p( 0 ) + v/i£^Jp( 0 ) + 


CM 


\p^d\ dr, 


(-4) 


ICol dr < e|(Co,Ci)(0)| + + Ve£A4Jp(0)) + l^ldr, 


whereas using allows to get directly 


P^ \d\dT <\{pb,d){0)\+ CM |ji|dr. 


Using the definition of ji and, again. Inequality p4.58p . we may replace ji with (Ci as follows 

(4.61) f \d\dT <Up{0) + VeCMJpiO) + CM f |Ci|fir+ f IpColfir. 

Jo Jo Jo 
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Then plugging ()4.59p and (|4.60l) in (14.611) and observing that e <C £A4, we get 

f\d\ dr < UpiO) + VICMJpiO) 

Jo 


+ 


pCM 


£|(Co) Ci)(0)l + (^p(O) + VeTA4j7p(0)) + pe J \d\dT 


Because we assumed that £ S> e, the last term may be absorbed by the l.h.s. Using in 
addition that e^A4 < 1, we end up with 

fi 


p^ / \d\dT< UpiQ) + ^eCMJp{^). 


Then resuming to (14.591) . (I4.60p . we obtain 


(4.62) 


' TlCol 

Jo 


dr + CJH \Ci\dT<Up{0). 


Obviously, this inequality implies that 


(4.63) 


CM / gi\dT<UpiO) + VICMJp{0)+e\(CoZm\- 


Of course, we get the same inequality if replacing and ("i with jo and ji. So one can 
conclude that for y/^Jn < p < cCM, we have (|4.58l) and 

(4.64) p^ f \d\dT + CM[ \ji\dT + pf |jo| dr < |(p6,d)(0)| + Ve£A4|(jo,4i)(0)|. 

Jo Jo Jo 

4.3. High frequencies. We eventually come to the proof of decay estimates for p > cCM, 
where c is some given positive constant. We shall use that fact that the systems satisfied by 
{b,d) and by {jo,ji), respectively, tend to be uncoupled for p —>• +oo. As regards (5, d), the 
classical approach for the barotropic Navier-Stokes equation, based on the study of 

Ul := 2|(6, d)|2 - 2pRe (hd) + \pb\^, 


guarantees, if p > c, that 

(4.65) |(p6,d)(t)|+p^ / |d|dT+ / p|6| dr < |(p6, d)(0)| 

do Jo 


Next, from the system fulfilled by (jo, ji), we get 


1 d 


2 dt 


idor + ijir +-ijor + 




+ 




CM 


n Jo 


\ji\ dr. 


d 


C. 


CM.^02 ic-\ 

|ji| =VrajRe( 6jo), 


, Re(jo ji) + -(1 + A4)Re(jo ji) + -^|jip-^|joP = \/n-Re(6ji). 

dt £ e-v/n £\/n e 


Therefore, for any k > 0 
1 d f ,2 1^ i2 ‘JkCM 

2s(w +l"l - — 




Re I 


/ . A 

.C.^ 

. .n CM 

( K \ 

1 + -^M 

- ?oF +- 


V > 

£ 

£ 

V y/n) 


kC‘^M 




(1 + A4)Re (jo ji) + Vn -Re 6 jo - 


- kCM 


P 


■Ji 


pe £ 

For p > cCM, it is clear that our choosing k small enough implies that 

• the first term of the r.h.s. may be absorbed by the second and third ones of the l.h.s.. 
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• we have |jol^ + liil^ - ^^^^Re(jo ji) « |ioP + 

. we have (l + + (1 - > f ^(UoP + liil"). 

Therefore, we end up with the following inequality 

(4.66) e\{jQ,ji){t)\+CM [ |(jo, ji)| dr < e|(jo,ii)(0)| +/I / H dr. 

do do 

Combining with (|4.65p . we conclude that if 

(4.67) p > cmax(l, £dV4) and p>CC, 
for a large enough constant C then 

(4.68) \{pb,d){t)\ + f \d\dT+ f p\b\dT <\{pb,d,£jo,eji){0)\, 

do do 

(4.69) |(io,di)(i)l + [ l(do,di)|dr < —|(p6,d)(0)| + |(j/o,ji)(0)|, 

e do pe 

whence 

(4.70) |(7'6,d,jo,di)(i)l + 7*^ / |d| dr + / p|5|drH-/ |(jo,di)l dr < |(p6, d, jo,di)(0)|. 

do do ^ do 

The only case where the condition p > CC may be stronger than p > cCM is when M 
is bounded. From our study for small p’s, we must assume that CM > and thus cCM 
is still much larger than \/2/n. Therefore, one may take advantage of (14.581) to bound the 
r.h.s. of (14.661) . and combining with (I4.65p . we get for p k, C 

(4.71) |(p6,d, \/e£jo,\/^i)(7)|+p^ / |dl dr + / p|6| dr 

do do 

+ ^ [ l(do,di)|dr < |(p6,d, \/eZjo, V^i)(0)|. 
do 

5. The non-equilibrium diffusion regime 

This section is devoted to the study of the so-called non-equilibrium diffusion asymptotics. 
Assuming that for some k > 1 and m > 0, we have 

(5.1) - > — and C^Cg —>• 

e nv 

we want to prove the convergence of the solutions of (jl.6p to those of (j2.5p or (j2.6p if m < -|-oo 
or m = -|-oo, respectively, when e goes to 0. 

The hrst subsection concerns the proof of global existence with ‘uniform’ estimates for the 
radiative Navier-Stokes equations (II.6h in the asymptotic (15.ip . in the small data case with 
critical regularity. Next, still for small critical data, we establish the global existence for the 
limit systems (|2.5h and (12.6p . In the last part of the present section, we combine the uniform 
estimates with compactness arguments in order to justify the convergence of (II.6p to (j2.5h 
or (12.6p . 
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5.1. Global existence and uniform estimates for (jl.6ll with £ « e and C?C.s > 1. In 
order to get a global-in-time existence statement for (11.611 in the non-equilibrium diffusion 
regime, we first put together the estimates that we obtained in the previous section, in 
the case £ ~ e and £^£s > 1. Even though localizing the linearized equations by means 
of Littlewood-Paley operators allows to get essentially optimal estimates for the linearized 
equations of (II.6p . it is not enough for our purpose, owing to the convection term u-Vb that 
may cause a loss of one derivative. The difficulty may be overcome by paralinearizing the 
whole system, as explained below. After that, it is easy to prove global in time estimates 
for the solutions to (11.611 just by combining the estimates for the paralinearized system, and 
standard product laws in Besov spaces to handle the other nonlinear terms. 

5.1.1. Linear estimates. Performing the change of variables (j4.1|l reduces the study to that 
of the linear system (14.4h . 

Low frequeneies estimates. Using (|4.28p . the comment that follows, (I4.39P and the fact that 
\{b,d,jo,ji)\ ^ |(b,h,jo,ji)|, we get for the solution {b,d,jo,ji) to (|4T]l, 

(5.2) \(b,d,jo,ji){t)\ + f \(b,d,jo,ji)\dT + f \]o\dT + M[ |ji|dr 

Jo Jo Jo 

< C\(b,d,jo,ji)iO)\ for all 0 < p < Ci, 
withQ To := jo - I and Ti := ji - ^^3o + 

Middle frequencies estimates. If £^£s ~ 1 then using (I4.54p . and the definition of Ci versus 
that of ji, we get 

(5.3) \{pb,d,jo,ji){t)\ + p'^ f \{d, jo)\dr + p f |6| dr-F A4 / |ji|dT 

Jo Jo Jo 

< C\{pb,d,jo,ji){0)\ for all Ci < p < c£A4. 

If C^Cg —)• - 1-00 then (14.4711 and (|4.48ll ensure that 

(5.4) \{pb,d,pjo,ji){t)\+p^ [ \d\dT + pf \(b,jo)\dr + Jii f |ji|dr 

Jo Jo Jo 

< C\{pb,d, pjo,ji){0)\ for all \/l + n~^ < p < er'VM, 

and, according to ()4.54ll 

(5.5) \{pb,d,{l + ^^^)jo,ji){t)\+ f \pb\dT + p^f \d\dT 

Jo Jo 

+ (^ + 7^)/ \pjo\dT + JA f |ji|dr < |(p6,d, (l-7^^)jo,ji)(0)l, 

Jo Jo 

whenever Cieuy/JA < p < cevfA. 


^Note that the last term of jo given by (14.221) is negligible for p £ 1 and may thus be omitted. 
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Hence 


(5.6) |(/) 6 ,d,max(l,min(p ,/9 ^e^uM))jo,ji){t)\ + p -^)jo)\dr 




+ p \d\dT + M \h\dT < C\{pb,d,max{l,mm{p,p e^i^M))jo, ji){0)\. 


High frequencies estimates. Using ()4.70p . we have 


(5.7) \{pb,d,jo,ji){t)\+ p^ 


f\d\ 

Jo 


dr + p | 6 |hr + A4 / \{jo,ji)\dT 


<C\{pb,d,jo,ji){0)\ for p>cCJA. 

For notational simplicity, we shall slightly abusively change c and Ci to 1 in all the 
following computations. 

Optimal estimates in Besov spaces. If C^Cg ~ 1 then localizing (14.411 with nonzero source 
terms / and g according to Littlewood-Paley operator A^, using (14.5p and following the 
computations leading to ()5.^ . ()5.3p and ()5.7p (combined with Fourier-Plancherel theorem) 
we end up for all s E M witlu 


(5.8) ||(u,jo,Ji)(t)b|,^ + \m\\%\^ + IIKOIlg+i + 


2,1 

ct 


111 rjs+2 d'T 
-^ 2,1 


+ 


/ ll( 5 ,jo,Ji)||Ji +2 fir + f || 6 ||^ 7 Vidr+ f ||jo|||i dr 

Jo ^2,1 Jo ^2,1 Jo ^2,1 


+ M 


h\fj!:^dT + 


'B. 


2,1 




B: 


< 


ii(s,jo,ji)(o)ii™+ii6(o)ii|; + ii(.(o)iiy*. + [' 

2,1 « 2,1 ^ 2,1 Jo 


e.i 


B: 


+ 


2,1 


h,l 

ds +1 

^ 2,1 


+ Il5b|jfir, 


with 

(5.9) ]o := jo - y/nh - y/nLdiY u and )i ■= ji +—Vjo - -^-3—V5. 

C yfJxCM CCsM 

According to our previous work in [7], the critical regularity framework corresponds to s = 
n/2 — 1. Therefore, the following quantities will play an important role 


\\{b,u,joJi)\\x ■■= + 

^ 2,1 


+ ll(^,jo,ji)IL5- 
^ 2,1 


and 


\\{b,u,jo,ji)\\Y ■■= sup \\{b,u,jo,ji){t)\\x + 

t>0 


hi 

n 

2 =’ 2,1 


+ 1 


+ 


+ 


ill 5_i + llloll + 

^2,1 ^2,1 


fef’l 

m,l,ei2A4 




dr 


..^+1 +MH(jo,ji)t:n''^)dT. 
^2,1 ^2,1 

We denote by X and Y the corresponding functional spaces (where time continuity is im¬ 
posed rather than just boundedness) and agree that Y(t) stands for the restriction of Y to 
the interval [ 0 ,t]. 


^Further explanations on the method will be supplied to the reader in the next paragraph. 
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In the case —>■ +oo, we have to change slightly the definition of the norms || • ||x and 

||y as the middle frequencies obey (15.6p . Consequently, we change || • Hx to || • ||xoo with 


\\ib,u,jo,ji)\\x^ ■= 


+ 


i£,l 


+ 




2,1 


n 

b7 


+ 


\h,l 

n 

^2,1 

2^5max(l, 

^2,1 l<2l^<eM 


2 '',min(2 VuAl))||Afcjo||L 2 


and II • ||y to 


\\{b, U,joJl)\\Y, 


sup \\{b,u, jo Ji)it)\\xg^ + 

t>0 



n 

^ 2,1 


dr 


+ 



,e,euM 


B. 


r-l 


+ 


2,1 



||^.’ 5_1 + M\\{joJi)t:^^^)dT 

Bl, / 

^ 2^t min(l,2^^e"^i/"^Al"^)||Auo||L2C?'^- 

1 < 2 '=<£AI 


To prove global estimates for the nonlinear system (II.6h . the natural next step would be to 
take advantage of (15.811 with s = n/2 — 1 and all the nonlinear terms in the r.h.s. Unfor¬ 
tunately, this does not work for the convection term u • V5 causes a loss of one derivative 

.71 .71 _^ 

(indeed, if b is in B^i then u ■ V6 cannot be smoother than B^i )■ A nowadays standard 
way to overcome the difficulty is to paralinearize (11.61) . that is to add to (14.41) the principal 
parts of the convection terms. This is the aim of the next paragraph. 


5.1.2. The paralinearized system. Before introducing the paralinearized system associated to 
(|1.6I) . let us shortly recall the definition of the paraproduct, according to the pioneering paper 
[2] by J.-M. Bony. The (homogeneous) paraproduct between two tempered distributions U 
and V satisfying (1,4. ip is given by 

TuV :=^Sk-iUAkV with Sk-i := 
k 


We also introduce 

T^U :='^Sk+2VAkU, 

k 


and, observe that, at least formally 


UV = TuV + T{^U. 


To some extent, if U is smooth enough then TjjV may be seen as the principal part of the 
product UV. This motivates our considering the following system 

dtb + ■ V6 -|- div u = F, 

dtu + Tff-Vu-Au + Vb- = G, 

(5.10) ^ 

^ divji -k f (jo - v^6) = 0, 

^ ^ Vjo + = 0 , 


where A = ;uA + (A-|-/x)Vdiv , M. := l + Cg, v and G are given time dependent vector-fields, 
and T is a given real valued function. 
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Proposition 5.1. For any smooth solution we have the following a priori esti¬ 

mate i/ 0 < m < +00 


\\{b,u,jo,ji)\\Y(t) < C*! ji)(0)||A + / \\Vv\\Loo\\{b,u,jo,ji)\\x dr 


+ 


f\\{VF,G)f%_,dT+ f\\{F 
Jo Jo 


-T^-Vb,G-T^-Vu)f\ dr 


^2,1 


A similar inequality holds if m = +oo, with Xoo{t) and Yoo{t) instead of X(t) and Y{t). 
Proof: Localizing System (IS.lOp by means of yields 
dtAkb + Ak{Tff ■ Vb) + dW AkU = AkF, 
dtAkU + Ak{T^ ■ Xu) - AAkU + XAkb - ^Akji = AkG, 
dtAkjo + ^ div Afcji + §{Akjo - VnAkb) = 0, 
dtAkji + ^ vAfcjo + = 0. 


(5.11) 


The important point is that in order to obtain all the estimates corresponding to p > Gi, 
one only has to resort to combinations between fluid unknowns on one side, and radiative 
unknowns, on the other side. This will enable us to use exactly the same energy method for 
(j5.10p as for (j4.4|) . in the middle and high frequency regimes, without introducing unwanted 
parts of convection terms in the inequalities. 

1. Low frequencies: 2^ < Gi. 

Including the para-convection terms in the source terms of (j5.1ip and repeating the compu¬ 
tations leading to (15.2h . we get after taking norms and using Fourier-Plancherel theorem 


(5.12) \\Ak{b,u,joJi)(t)\\L2 + 2' 


.2k 


f 


\Ak{b,u)\\L2dT+ ||Afcjo|A2dr 


+ k'Cs 


WAkhh^ dr < ||Afc(6,u, jo,ii)(0)||i,2 


+ ^ \\Ak{F-T^-Xb)\\L2dT+ I \\AkiG-T^-Xu)\\L2dT. 


2. Medium frequencies: Gi<2^ < cCCg- 

Keeping in mind the proof of (j5.8h . we see that it is suitable to introduce 


Cl := ji + 


1 


-Vjo- 


y/nCM 

Now, because we have 

J dtAkb + div AkU = AkF, 

1 dfAkU — AAkU + vAfcft = ^^vAfcjo -|- ^^AkCi + AkG, 

we easily get by computing 

(5.13) — — ^2||(Afc6, Afcu)||^2 + ||AfcV6||^2 + 2(AfcV6| Afc'u)2,2^, 
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and by using Lemma 4.1 in to handle the para-convection terms, the following inequality 
for all 2^ > Cl 


||(AfcV 6 , Aku){t)\\i2 + 2 


2k 


Jo 


L2dT+ / ||AfcV 6 || 2,2 dr < ||(AfcV 6 , Afcu)( 0 )||i 2 




f ||(AfcVF, AfcG)||^2 dr-F f llvAfcjo 11^2 dr£A 4 [ || VA^ji 11^2 dr 
Jo Jo Jo 




iVdllroo II(AA:/V 6 ,Afc/'u) 112,2 dr. 


Then looking at the equations satisfied by A^jo and A^i^i (in the spirit of (|4.40p L we derive 
inequalities similar to (I4.50p and p4.5ip for ||Afcjo|| 2,2 and ||Afc(Ci|| 2 , 2 , and thus following the 
computations leading to p4.54l) . we end up in the case m < -|-oo with 

(5.14) ||Afc(V6,M, Vjo,Ji)(t)||L2-h 2^^ / ||(Afcn,Afcjo)||L2^^T+ / ||AfcV6||2,2 dr 

Jo Jo 

l|AfcCi||L2dr<||Afc(V6,d,VjoJi)(0)||i2+ / ||Afc(VF,G)||^2dr 

Jo 

+ Y. [ l|V^T||Loc||Afc,(V6, 


+ k/Cg 


u)\\l2 dr. 


Comparing the definition of Ci and ji, we see that one may replace Ci with ji above, if 
Cl < 2*^ < cCCg. 

The obvious modification^ to be done if m = -|-oo are left to the reader. 

3. High frequencies: 2 ^ > cCCg- 

Again, we compute (I5.18|l to bound the fluid unknowns. In addition, to handle radiative 
unknowns, we compute for some small enough k (see the proof of (I4.66|l l the following 
quantity 

“(l|Afcjo|li2 + ||AfcJi||i2 - kAA 4 2"2''(Afcjo|AfcdivJi)2,2). 

Combining the computations leading to (15.711 with Fourier-Plancherel theorem and Lemma 
4.1 in [7] eventually yields 

||Afc(V 6 ,d,jo,ii)(i)||L 2 -h 2 ^*^ [ ||Afcu|| 2,2 dr-k 2 ^ f \\Akb\\L 2 dr+uCg f ||(Afcjo, Afcji)|| 2,2 dr 

do do do 

ifc(V 6 ,n, jo,ji)( 0 )||L 2 -k / ||Afc(VF, G) 112^2 dr-k ^ [ ||Vi;||i,<=° ||Afc/(V 6 , n)|| 2,2 dr. 

Finally, multiplying (j5.12p . p5.14p and the above inequality by 2 ^^ 2 'and summing up over 
k completes the proof of the proposition. □ 

5.1.3. A global existence result. According to the computations of the previous paragraph and 
to the change of variables (14.ID . it is suitable to introduce the following norms for getting 
global solutions with uniform estimates in the cas^ m < -koo 

\\{b,u,joJi)\\x- ■■= \\b\f:''n_^+u\\bf:'i '-k ||(n, jo,di)ILf-i, 

"i? ^ R ^ Brf 


< IIA 


^Writing out the corresponding definition if m = +oo is left to the reader. 
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and 


ji)||y/ := sup \\{b,u, jo, ji){t)\\x’' + 
t>o 




B. 


Hi 


+ 


+ 


Bli 


)o|| . «-i + 
b7, 


\h,u 

n 

2,1 -^^ 2,1 
-1 


+ ^\\u\\.^+i 

^ 2,1 


dr 


m,u ^ 




B, 


Hi 


h,EM 


B^ 


with jo := jo-b - fdivn and ji := ji + -^^jo - 

Of course, if {b',v!, J'q, j'l) and {b,u,jo,ji) are interrelated through (14.ip and vC is used 
for {b',u', J'q, j[) instead of C, then we have 

\\{b',i^,joJ[)\\x} = i'~^\\{b,u,jo,ji)\\x- and ||(^', "S', ji)||yi = z^"^||(5, n, joJi)||y/- 


Theorem 5.1. Assume that T ~ 1, that liminf e“^nz^T > 1 and that C^Cg ~ 1- There 
exists a positive constant p depending only on fi/u, n and on the pressure law such that if 
e is small enough and the data ^oOo,oO'i,o) satisfy 

(5.15) m,ro,jl^oJlo)\\x^ <W, 

then System (HSl) admits a unique global solution (b^,u^, Jq, jf) in Tf. In addition, we have 

(5.16) 11(6^,HoH'i)l|y/ < C'||( 6 o,izo,jo, 0 )Ji,o)llA:|'. 

A similar result holds true if C^Cg H +oo. 


Proof: Performing the change of variables proposed in ()4.ip reduces the proof to the case 
z^ = 1 (changing C into C := uC). Hence we consider a smooth enough solution to ()4.2p . 
and show that one may close the estimates globalljH under Assumption ()5.15l) . 

Let us set Uq := ||(^o> ^>io,o>ii,o)llAi and U{t) := ||( 6 ^ jg, jf)||y^i(t). In what follows, 
we drop exponents £ for notational simplicity. Finally, to shorten the presentation, we just 
treat the case where C^Cg « 1 . 

Now applying Proposition 15.11 with v = u, F := —• u — ki{b)divu and 

CM 

G := ■u + k2ib)Au-k3ib)Vb+^ k^{b)ji, 

n 

yields for all t > 0 


(5.17) Uit)<C(Uo+ I \\Vu\\L^\\ib,u,jo,ji)\\x}dT 

rt 


'0 


+ 


'0 


|F|| . + ||F - ■ V 6 || + ||G|| m • Vn|| 5-1 ) dr 

-^ 2,1 ^ 2,1 -^ 2,1 -^ 2,1 


Using standard continuity results for the paraproduct and remainder, and composition esti¬ 
mates leads to 

\Ty^-u\\^n < C'||V6||^n_i||M||^-+i, 


'B. 


2,1 


'B. 


2,1 


ci( 6 )divn|| ™ < C\\b\\ ^ ||divn|| ™ . 

^ 2,1 ^ 2,1 -^ 2,1 


^Existence follows from spectral truncation as in e.g. [T], Chap. 10, and is thus omitted. As for uniqueness, 
we refer to [7]. 
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Hence we have 
(5.18) 

We also have 


/ 


IFII .n dr < CU^it). 

I " D ^ ' 

^2,1 


\Tit ■ < C'||n||^--i||Vn||^n , 


2,1 


'B. 


2,1 


'B. 


2,1 


IT' 


Vn • ^ C'llVdII II^IL^-O 

^2,1 ^2,1 ^2,1 


\\k2{b)Au\\ ^.r<c\\b\\ ^ liv'nil 

-^2 1 '^2 1 ^2 1 

||fc3(6)V6|| »_i<C||6|| . ||V6|Ln 


Bounding /lAl/i; 4 ( 6 )ji is slightly more involved as it is not true that the low frequencies of 


ji are bounded in L^(R^; B 21 )• However, one may write that 
£Mji = 

Therefore 
£A1||A:4(6 );i| 








LliBliY ^riB^i) 


+ \\b\\ 


L^Bli) 




'BUB?.i Y 


L^Bii y 


Hence 


(5.19) 

Allr,^ • Vn|| + ||G|| dr < CU\t) 

JQ ^ 2,1 -° 2,1 

Finally 

||d-V 6 || <C\\u\\ n ||V 6 || 

^ 2,1 ^ 2,1 ^ 2,1 

||fei( 6 )divn|| <C|| 6 || » ||divu|| n_i. 

-^ 2,1 -^ 2,1 ^ 2,1 


Therefore, by Cauchy-Schwarz inequality 


(5.20) 


f\\F-T^-Vb\\^.^dT<CU\t). 
Jo '°2,l 


Inserting (|5.18p . (I5.19p . (15.201) in p5.17p and remembering that B 21 '—>• L°° (to ensure that, 

say, | 6 | < 1/2 if ||5|| . ^ is small enough), we end up with 
B2,1 

U{t) < C{Uo + U^{t)) for all t > 0. 

By a standard bootstrap argument, we easily deduce that 

U{t) < 2CUo for all t > 0, 

provided the data have been chosen so that 4(7 ^Hq < 1 - C 
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5.2. Study of the limit system. In this paragraph, we prove the existence and uniqueness 
of strong (small) solutions with critical regularity for Systems (12.5p and (12.6p . We shall give 
a common proof that works for both systems. 

Before giving the global existence statement, let us introduce the solution space 

• If m G (0, +oo) (that is for System (12.5M then Initial data will be taken in the space 
which is the set of triplets {b,u,jQ) satisfying 


\\ib,u,jo)\\x-' ■■= 


" dT- 
^2,1 


+ ^l 


\h,u ■ 

n 

^2,1 


+ ll'ull . 'i-i + u 


^2,1 


£,u-^ 

n 

^2,1 




-1 


n o CX), 
^2,1 


and the solution space will be the set of triplets (b, u,jo) in Cb(M+; satisfying 


\\ib,u,jo)\\y‘'-=sup\\{b,u,jo){t)\\x-+ {\\jo-b\\ n +v 


t>o 




h,i^ 

n 


dr. 


If m = +00 (that is for System (12.61) 1. Initial data will be taken in the space X^ 
which is the set of triplets (6, u, jo) satisfying 


V 

OD 


\\{b,u,jQ)\\x^ := 


^2,1 


^ 2.1 ^2 1 


2,1 


and the solution space y^ will be the set of triplets (b, u, jo) in Cfe(M+; X^) satisfying 


\\{b, u,jo)\\yi^ 


sup||(5,u,jo)(t)||;tii 

t>0 


+ 







+1 


+ 


h,v 

n 

^2,1 


dr < oo. 


Theorem 5.2. There exist two positive constants c and C so that if 


(5.21) ||(^o,iio,^0,0)11^- < CB {case m <+oo), 

(5.22) or ||(^o,ifo, jo,o)||A^ < CB (case m = +oo), 

then System (|2.5I) (resp. (12.61) ) admits a unique solution in the space y^ (resp. y’fa) satis¬ 
fying in addition, 


(5.23) ||(6,tt, jo)||y- < C'||(&o,no, jo,o)||A‘' if m <+oo, 

(5.24) II(5, u, jo)< <711(50,^0, jo,o)II if m =+oo- 

Proof: Set k := k/h and fh := mn. As usual, it is enough to treat the case b = 1 as 
performing the change of unknowns 

{b,u,jo){t,x) = {b,u,jo){n~H,u~^x), 

gives Systems (|2.5I) or (|2.6p for (6, u, jo) with b = 1 and A := Ajv and, obviously 
||(6,u, jo)(t)b- = n\\(b,u,jo){iy-^t)\\xi and ||(6,u,jo)b- = ^\\(b:U,jo)\\yi. 

Let us start with the study of the linearized equations with no source term, namely 


(5.25) 


dtb + div n = 0, 
dtu — Au + V6 + n“^Vjo = 0, 

, dm + K(jo - Ajo - 6) = 0. 
The divergence-free part Vu of the velocity satisfies 


dtPu — p,AVu = (), 



DIFFUSIVE LIMITS FOR A BAROTROPIC MODEL OE RADIATIVE FLOW 


31 


(5.26) 


while the coupling between b, d ;= A^^divu and jo is described by 

dtb + Ad = 0, 

dtd — Ad — Kb — n~^Aj() = 0, 

, dtjo + K{jo - fh-^Ajo -b) = 0. 

Note that the stability of a similar system has already been established in the previous section 
for K > 1 (or, equivalently, k > 1/re). 

Linear estimates for low frequeneies. We introduce Co := jo ~ ^ ~ K~^Ad and notice that 
dtb + pd = 0, 

dtd + ^)d- (l + i)p5= ipCo, 

<9tCo + («^ + (t + ^)^^)Co = -((1 + + §)p‘^b+ ((1 - 

On one hand, because Hn > 1, the method described in the appendix (see in particular (IB.71) 1 
allows to write that, omitting the dependence with respect to k 

\(b,d){t)\+p^ f \(b,d)\dT <\(b,d){0)\+p f \Co\dT. 

Jo Jo 

On the other hand, the last equation directly gives 

\Coit)\ +(k + + ^)p‘^) [ \Co\dT <\CoiO)\ + c(l +—)p‘^ [ \(b,pd)\dT. 

\ m Kn / Jq ^ rer/ jo 

Hence plugging the second inequality in the first one 

i(6,d)(t)iv f\{id)\dT<\{id)m+ 

Jo 


p 


l + (l+m 


Tj^(^IO)(0)| + (l-Kre ^)P^ \(b,pd)\dT 


l+m 


Hence 


It is clear that the last term may be absorbed by the integral of the l.h.s. if p <C 
we eventually get for some small enough > 0 

(5.27) |(6,d,pCo)(t)|+p^ / \(b,i)\dT+f |pCo| < |(6, d, pCo)(0)| if p < ( ^ )pe- 

Jo Jo Vl + rre/ 

Linear estimates for high frequencies. We set 6 := d — A“^6 and notice that 

dtb + b = -p5, 

dtd + (p^ - l)d = p~^b + re"V/o, 


Therefore 


. dtjo + k( 1 + jo = Kb. 

I^(^)l + {P^ “ 1) / 1^1 ^ l^(0)l + P~^ [ \b\dT + - f p|jo| dr. 

Jo Jo n Jo 


At the same time 


\m\+ t\b\ 

Jo 


dr < |6(0)| + p / |d| dr. 


|jo(t)| + re( 1 + —) / |jo|dr < |jo(0)| + re 

Jo 


f\h\ 

Jo 


dr. 
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Hence 


\m\ + iP^- 


ip^-l) [ |?|dr<|?(0)| 

Jo 


+ (i + -)|6(0)l 

p n/ 


n ^ 


+ ( )pb’o(0)| + ( H-/, I ~-l 2', 

Vl+m ^/9^/ V n(l + m ‘^p^) 


|5| dr. 


Therefore there exists a constant ph depending only on m and n (with n > 2 if m = Too) 
such that for p > ph, we have 

(5.28) \{pb,6){t)\+mm{p,mp~^)\jo{t)\ + p [ \(b, jo)\dr + p^ f \d\ dr <\{pb,6, pjo){0)\. 

Jo Jo 

Of course, one may replace 5 with d in (|5.28l) . 


Linear estimates for medium frequencies. The stability argument used just below ()4.37p 
allows to write that there exist two constants c and C depending continuously on 1/m, such 
that if p G Ph] then 


(5.29) 


\{b,d,jo){t)\ <Ce |(6,d, jo)(0)|. 


Estimates for the paralinearized system. The previous steps allow us to get handy estimates 
for the following paralinearized version of System (12.61) 


(5.30) 


dtb + Tff ■ Vb + div u = F, 
dtu + • Vu — Au + V6 + Vjo = G, 

. dtjo + K{jo - m~^Ajo -b) = 0. 


More precisely, following the steps leading to (15.271) . (j5.28p and (|5.29p . introducing (q := 
jo — b — K“^divu, and arguing as in Subsection 15. 1.21 we end up wittH 

(5.31) \\jo{t)\f:\ +\\{b,u){t)\f:\_^+ j [\\jo-b\{-^\ +\\{b,u)\{-'\ )dT 

R^ In R^ R^ 

^ 2,1 ^ 2,1 U ^ 2,1 ^ 2,1 


< 

r\-i 


IlioWir:', +ll(6,a)(0)ir:V, + 


e,i 


^2,1 


^2,1 


|T - • Vb\ff\_, + \\G - • VulffiAdr. 


'0 


B^~ 

^2,1 


B^' 

^2.1 


For high frequencies, we get, in the case m = +oo 

(5.32) \\{b,jo)it)t:\ +\\uit)f:\_^+ [ {\\{b,jo)t:\ +\\ut:\+,)dT 

B^ B^, Jo B^, Bl^ 


< 


ll('>,io)(o)f;l, + ||s(0)rr.. + 


h,l 


B 2 
-^2,1 


B 2 
-^2,1 


r% +\\Gt.%Adr 
^2,1 -^2,1 


+ / l|Vi;||i,oo(||(6, jo)|| + ||u|| n_i) dr, 

Jo ^2,1 ^2,1 


5 


Here we do not track tire dependency with respect to m. 
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and if 0 < m < +oo 

(5.33) \\b{t)f:i + lbo(i)||^’n_2 + + [ {\\{b,jo)\\% +\\u\\%+i)dT 




^2,1 


^2,1 


10 


< 

rs_/ 


iii>(o)ii'‘4 + 

^2,1 


h,l 

n 

^2,1 


+ ll^(0)ll^.’5-i + 

^2,1 


^2,1 


\Ff% +\\G\H,]dT 
^2,1 ^2,1 


^2,1 


+ 


[ l|Vn||L-(||6|| 

Jo 


So 


+ 


h 

n 

B^~ 

^2,1 


+ \\u\\^^-i) dr. 


2,1 ^ 2,1 2,1 

Proof of existence. We only establish global-in-time a priori bounds in the space 3^^ or 3^^^ 
for the solutions to (j2.5p or (j2.6p with data satisfying (I5.2ip or p5.22p . Our proof is based on 
P5.3ip . (j5.32p and ^5.331) with v = u 

F = —Tyj, • i; —/ci(6)div ft and G =■ u + k 2 {b)Au — ks{b)Vb — n~^k 4 ^{b)Vjo. 


Bounding \\F — Tjj ■ Vb\f ]\and ||F||^’n relies on (j5.18p and ()5.2np . As regards G, the 


I 


So 


Bo 


computations that we did in the proof of Theorem 15.11 ensure that the hrst three terms may 
be bounded as in (j5.19p . To handle the last term, k 4 ,{b)Vjo, in the cas^ m < -|-oo we use 
the decomposition 

h{b)Vjo = ki{b)Vb + ki{b)V{jo - b). 

The hrst term may be bounded quadratically exactly as ko{b)'Vb. As for the last term, we 
may write 

.f J|V(jo - ^)||, 


||A:4(6)V(jo-6)|L,.t < 

v-^2,1/ 




"L-(B2^i)' 

hence it is also bounded by G\\{b,u, jo)\\yi(^ty 

This enables to conclude that we do have for all t G M_|_ 

\\{b,u,jo)\\y^t) < C{\\{b,u,jo){0)\\x^ + ||(&,w,Jo)||yi(t))- 
This obviously yields (15.241) if (15.221) is fulhlled. 


Proof of uniqueness. It works the same as for the standard barotropic Navier-Stokes equa¬ 
tions: we look at the system satished by the difference {S),5u,8jo) between two solutions 
{b^^u^Ao) [b^^d^Go) of (I2.5h . namely (denoting ATj = 1 -|- A:, for i = 1,2,3,4) 

' dtdt) + ?22 . V& = -Su-Vb^ + {Ki{b^) - A:i(62))divM2 _ Ki{h^)6:iYSu, 

dtSu + u^-VSu + Su-ViA - {K2{b‘^)-K2{b^))AiP - K2{h^)ASa + {Ko{b‘^)-Ko{h^))Vb‘^ 

' +A3(61)V* + n-\K^{h^)-Ki{b^))Vjl + n-^K^QA)V5jo = 0, 

, dtbjo + k((5jo - & - iAhjo) = 0. 

Now, exactly as for the barotropic Navier-Stokes equations, it is possible to bound &, Su 
and djo just resorting to basic estimates for the transport and heat equations. However, 
the hyperbolic nature of the hrst equation forces us to estimate {S),du,^o) with one less 
derivative, namely in 


L°°(0, 


• 2--1 


n _ ry n_ 

(L-(0, T; % ) n L\0,T-, x L\0, 


. IL_1 

T-.Bfy ). 


In dimension n = 3 combining estimates for the transport and the heat equation allows to 
get uniqueness on a small time interval, then on the whole M+ by induction. In dimension 


®The case m = -l-oo does not require that decomposition. 
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n = 2 , this is slightly more involved as some product laws do not work correctly if estimating 
{6b, 6u, 6jo) in the above space (some regularity exponents become too negative). Nevertheless 
this may be overcome by combining logarithmic interpolation and Osgood lemma (see e.g. 
[ 6 ] for more details). This completes the proof of the theorem. □ 

n 

Remark 5.1. If 0 < m < +oo then one may alternately assume that jo is in Taking 

advantage of the parabolic smoothing given by the equation for jq, it is not difficult to get a 
solution ( 6 , 12 , jo) with 



5.3. Weak convergence. Here we justify the weak convergence of (II. 6 p to (j2.5p or (12. 6 p 
under the assumption that liminfe^Ts > 0 and that £ tends to ^ for some k> 1. 

Theorem 5.3. Let the family of data ( 69 ,1^, jo oO'i o)o<£<i satisfy Condition (15.151) . As¬ 
sume in addition that 

(5.34) —)• m G (0,+ 00 ] and ^ a € (1, + 00 ). 

Then the global solution {b^,u^, Jq, jf) given by Theorem \5. 1\ satisfies 

, n_-| . n 

jf —)■ 0 in Bf^^), 

and, if {bff,uff, Jq q) (bo,uo,Jo,o) then {¥,vF,jQ) converges weakly to the unique solution 
( 6 ,12, jo) of (|2.5I) supplemented with initial data ( 60 , l2o, jo,o) • 

Proof: From (j5.16p we gather that 


= 0{M-^) and = 0{M-^) in L^(R+;H|i ). 

Therefore, taking advantage of the boundedness of the low frequencies of Vb^ and Vjg in 

n 

Li 

we get 
(5.35) 


Bfi ), and of the fact that 


fi = fi 


1 


CM 


V jo" + 


1 


CCsM 


--1 


V 6 ^ 


jl = OiE) in L\M+,Bl, +BI,). 


^ . IL_1 

Next, we observe that (15.161) implies that {¥) and (u^) are bounded in L°°{M.^; Bf-^^ (H 

. IL_|_1 . H • IL _1 . 

Bfi) Cl L^{M.^; Bfi + Bfi) and L°°(K^-,B 2 i ) H L^(M_|_; ), respectively. Note that 

• IL 

this implies that vC is bounded in L'^{]^^; Bfi)- Because 

dtb^ = -if ■ V6^ - A;l(6^)div^2^ 


, n_ 2 . IL • IL _ 

and the product maps Bf^ x Bf^ in Bf^ , we thus get in addition that dtV is bounded in 
Bf I ), and thus {¥) is bounded in {M._^_■, Bfi ). Interpolating with the bound in 

■ C a ■ C—a 

Cf,(M+; Bfi)^ (^’") bounded in C 2 (M+; Bfi ) for all a G [0,1]. Then combining 

locally compact Besov embeddings and Ascoli theorem allows to conclude that there exists 6 
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in B 21 ^ 21 ) ^ ^21 + ^ 21 ) ^ sequence (efc)fc6N going to 0 so that, 

for all 0 G 5 and all a G (0,1] 

n _ 

(5.36) in LToM+^BI^ )• 

From (jS.lGD . we readily get for some sequence {ek)k&N tending to 0 

(5.37) in L°°(M+;B|i )nLi(M+;5|^ ) weak*, 

which, combined with (|5.36|) is clearly enough to pass to the limit in the mass equation. 

. 21 — 1 

Next, we see that (|5.16p implies that (jq) is bounded in (R+; ). Hence there exists 

• -—1 

jo G L°°(M+; ) and a sequence {ek)keN going to 0 so that 


Because 

(5.38) 


■ €]c 

Jo 


Jo 


in L° 




■ 21 — 1 

B 21 ) weak *. 


1 


-div jf =- 


B f 




^fo - 


C{\ + Cs 


-5tdivjf, 


and (|5.35l) implies that dtjl —)• 0 in the sense of distributions, we deduce that 

1 , A W 

-divjf -Ajo m S . 

£ nm 

Note that the right-hand side is 0 if m = -|-oo. Therefore (6, jo) satisfies the third line of 
(j2.5p (case m < -|-oo) or ()2.6p (case m = -|-oo). 

Let us finally pass to the limit in the second equation of (II.6p . The main difficulty is 
that, owing to the radiative term which is only bounded in a -in-time type space (namely 

L^(M+;i72^^ ) or so), one cannot take advantage of some suitable bound of dtiF so as to 
glean some equicontinuity and then resort to Ascoli theorem. To overcome this, we use the 
fact that, owing to p5.38p 


dt -h ^(1 = -if ■ Vtf -h (1 k2{V))Aif 

-(1 + kom)VW - - (1 + ki{W)) Vjo" + -k'^{W)dtW]l. 
n n 

,77. _^ ,77 

Now, because {if) is bounded in the space (M+; ^ 2 ^ ) n L^(M_|_; 7^2^^), (jg) is bounded 

,77_-| . n 

in L°°(R^; B 21 ) (b^) is bounded in (L^ n L°°)(M+; H 2 ^^), product laws in Besov 

spaces ensure that the first four terms of the r.h.s. are bounded in L^(R+; B 21 ) (only 
in L^(R^; B 200 ) if R = 2). The same property holds true for the last term for {dtV) is 

bounded in LfiRj^-, B 2 i ) and (jf) is bounded in L°°(M_|_; 732^]^ ). Using locally compact 

Besov embedding and Ascoli theorem, one can now conclude that there exists some v in 
■ --1 

L°°(M+; B 21 ) so that for all 0 in 5 and a G (0,1), we have, up to extraction 
+ ^{'^ + ki{b^))j{j — >(j)v in Lf^^{R+; Bi~^~°‘). 

Of course, combining with p5.35p . we discover that v = u. Hence we also have 
(jif —> 4>v in Lf^^{R^B 2 ^i ) for all 4> ^ S. 
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It is now easy to conclude that the second line of (12.51) is fulfilled by (5, ft, jo). 

Finally, that the whole family ( 6 ^, ft^, Jq) (and not only subsequences) converges to {b, u,jo) 
stems from the fact that the solution to (I2.5p or (12. 6 p is unique. □ 

Remark 5.2. It is also possible to justify the strong convergenee of the solutions of ()1.6p to 
(ESP or dMP using (|5.35p and performing the difference between ( 6 ^,'!F,Jq) and the solution 
( 6 , ft, Jo) to the limit system. Again, taking advantage of the deeay properties o/jf is crueial. 
Note however that, exactly as in the proof of uniqueness, owing to the hyperbolic nature of the 
density equation, one eannot prove the strong eonvergenee in the solution space. There is a 
loss of one derivative that may he partially compensated by combining with uniform estimates. 
As we do not think this approach to bring much compared to weak compactness, we leave the 
details to the reader. 


6. The equilibrium diffusion regime 

This section is devoted to the mathematical justification of the equilibrium diffusion regime 
given by ()2.inp . To avoid useless technicality, we focus on the case where 

( 6 . 1 ) C —)• +00 and eCM. ~ 1. 

6.1. Linear estimates. Let us gather the estimates we proved for ()4.7[) for the above asymp¬ 
totics in Section m 

Regarding low frequencies, one may combine (j4.20p and (14.2111 to get 

___ ft ^ ^ 7 ft C.AA r 

(6.2) \{h,d,jo,ji){f)\+p^l \{b,d)\dT + - I \}o\dTl -/ |ji|dr 

Jo ^ Jo ^ Jo 

<C\(b,d,jo,ji){0)\ for 0 < p < \/l + n-^, 
with Z := vL, jo := jo - 'jZb - ^ ^pd and ji = ji - ^^6. 

For middle frequencies, we have according to ()4.58l) and ()4.64l) 

(6.3) \{pb,d,jo,ji){f)\+ [ \pb\dT + p^l \d\dT + p ( |jo|dr + £A4/' |ji|dr 

Jo Jo Jo Jo 

< C'KA^ jo, ji)( 0 )| for y/2jn < p <cZjA, 

and (14.701) gives, if A4 is large enough 

__ ft ^ ft ^ 7 ka ft ^ ^ 

(6.4) \{ph,d,jo,ji){f)\+p^l \d\dT + pl \b\dT^ -/ |(jo,ji)|dr 

Jo Jo ^ Jo 

<C\{pb,d,jQ,ji)Zf)\ for p>cZjA. 

If A4 is bounded then we must assume that p > CiCAi for some C\ > c. However, we have 
(j4.7ip and A4 bounded implies that e£ « 1. Therefore (16.3p is satisfied up to p < CiCJ^. 
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For the whole system (j4.4p with nonzero source terms / and g, we thus obtain (taking 
slightly abusively c = Ci = 1 for notational simplicity) 


(6.5) ii(ajo,ji)(*)iiB,+iii>(«)iri; +iiK*)ii^.ti +/ 

2.1 1^2.1 1^2,1 Jq 

[ llhlll! dT+ f \\{b,joJi)\f^]+2dT + j IIj 

Jo "2.1 Jq "2.1 Jq 


u\\^s +2 dr + 


C 


2,1 


11^.1 ^ 
o|L. dr 


CM 


^I+I dr + CM 

"2,1 


£ Jq " ' "2,1 

■t 


/ 


Ji\\c. dr 


'B. 


2,1 


+ 


/Vi4’.Vi 

Jo ^2,1 


dr 


CM 


\\{jo,3i)tj!:^ dr < ||(n,jo Ji)( 0 )||, 


I JDS 
^ 2,1 


2,1 


+ 


W0)l|U^+l|6(0)||^..V+ r(||/||6‘_ 


+ ll/ll^’s+1 + 

^ 2,1 


\B. 


(.) 


with 


]o'■= jo — Vnb — ^/n ^d\v u and ji = ji + V 6 . 

C CM 


Back to the original variables, that linear analysis induces us to introduce the following norms 


ll(^,«,io,Ji)|| v. := 


' + \\u\\.^}-i + ||(io,Ji)|L^-i and 

"21 "21 "2,1 "2.1 


,u,joJi)\\y. ■■= sup ||( 6 ,m, jo,Ji)(i)|| V. + V 


i >0 

+ 


I - 
n 

"2,1 


+ 


/ 

4R- 


CM,,-, 
+ ^llh 

m,v~^ ,CM 

n 

1 ^ 2,1 


[ (\\{b,Jo 
Jm+\ 

\ 11 ^,^ 1 \\ -^\ 

,Jl)|| . -h \\U\ 

"2,1 

^ 1 

-t 

llioll^.% dr 

"2(1 e 


+ CM\\ji\j 

,CM 1 CJA 

n 1 ~r 

R 2 E 

^2,1 ^ 


^2,1 


+ 1 


dr 


,h,CM 

1 ^ 2,1 


with jo := jo-b - ^div-u and ji := ji + £^V 6 . 

We denote by and the corresponding functional spaces (where time continuity is 
imposed rather than just boundedness). Of course, we still have 

\\{b',ijf,jo,j[)\\xi = i^~^\\ib,u,jo,ji)\\x. and ||( 6 ', i/, Jq, Ji)||^i = V^||( 6 , n, jo,ii)||y., 


through the change of variables (j4.ip . if we replace £ by £ in the left-hand side. 


6.2. The paralinearized equations. In the equilibrium diffusion limit case the estimates 
for the paralinearized system 

dtb + ■ V 6 -|- div u = F, 

dtu + ■ Vu - Au + Vb- =G, 

( 6 . 6 ) { ^ 

. dji + ^ + = 0. 


recast as follows 
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Proposition 6.1. For any smooth solution we have the following a priori esti¬ 

mate for (j5.10p 

jo, Ji)lly^.(i) < c(^\\{b,u,jo,Ji)iO)\\xi + ||Vi;||Loo||( 6 ,'u,jo, ji)||jf. dr 

+ r ||(VF,G)f.’i_, dr+ r ||(F-r^-V 5 ,G-r^-Vu)||^:^„_, dT+ f \\T^-Vb-Frf^^^ dr) • 

do ^24 do do Bar ' 

Proof: Except in the middle frequencies range, the proof goes along the lines of the 
corresponding result in the non-equilibrium case. Still assuming that v = \ and replacing 
C with C = vC then, working directly on the localized paralinearized system (16.6p . and 
combining Inequalities (j6.2p to (16.4p with estimates for the para-convection terms gives 

1. Low frequencies: 2^ < Ci. 


fb,u,joJi){t)\\L2+2“^’^ \\Akib,u,joJi)\\L2dT + - / ||Afcji||L2dr 

do s do 


+- f ll^fcjo||L2 dr < ||Aa.( 6 ,m, do, ji)(0)||L2 


+ / ||Afc(P-rr-V 6 )|A 2 dr + 


fw^kiG 

do 


- T^-Vu)\\l 2 dr. 


2. Medium frequencies: Ci < 2^ < cCM.. 

One has to keep in mind that in order to derive (16.3p from ()4.58p and (j4.64p . one has to 
consider the system that is fulfilled by ( 6 , ft, Co, ji) with Co := jo ~ y/nb. In particular, a part 
of the the paraconvection term of b enters in the equation for Co as we have 

£ 1 ^ 

dtCo + —Co H-^divji - v^divM = y/n{T^ -Vb- F). 

£ £\/n 


Therefore, following the computations leading to (I4.58p and (I4.64h . and using Lemma 4.1 in 
[7] to bound the convection terms coming from the equations for b and u, we end up with 

||Afc(V6,u,do, ji)(i)||L22^'= / WAkuWL^dr+ 2^ f ||(Afc6, Afcjo)||L2 dr 

do do 

+C{l + Cs) [ ||Afedi||L2 dr < ||Afc(V6,u,do, ji)(0)||L2 f ||Afc(VF,(5)11^2 dr 

do do 

+ [ ||AArr-V6-F)|A2dr+V / ||VT||r-||Afc,(V6,u)|A2dr. 

3. High frequencies: 2^ > cCAi. We get 

||Afc(V6,u,do, ji)(i)||L2-h f {2^’^Aku\\L2-^2^\\Akb\\L2) dr + CMe f ||Afc(do, Ji)||i,2 dr 

do do 

Afc(V6,£,do,Ji)(0)||r2+ [ \\Ak{VF,G)h2dr 

do 

+ J2 [ l|VT||r^||Afc,(V5,£)|A2dr. 


< llAi 




Putting together all those inequalities completes the proof. 


□ 
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6.3. A global existence result. Our global existence result with uniform estimates reads 
Theorem 6.1. There exists a positive constant rj depending only on n and on the 

pressure law such that if £ & (0,1) and the data (^O’^oOo oOf o) satisfy 

(6.7) \\{bo,UoJoflJlo)\\x^^ < W, 

1 1 

then System (HSl) admits a unique global solution (b^^, u^^, Jq , jf) in Yf'. In addition, we have 

(6.8) \\{h^,If,jl,3l)\\y, < C||(6S,nS,jo,oo!,o)llx.- 

The proof relies on Proposition 16.11 Note in particular that the ‘new’ last term in the 
estimate of () 6 .ip does not entail a loss of derivative as we simply have 


r^-< llu-V 6 || , 


< 


The rest of the proof works exactly the same as in the non-equilibrium case. 


□ 


6.4. Weak convergence. Here we justify weak convergence to (j2.inp when assumption (j6.1|) 
is fulfilled. 

Theorem 6.2. Let the family of data o)o<£<i satisfy (16.7p . Then the global 

solution in Yf' given by Theorem \6.1\ satisfies 

Ti_-| n 

in +4p), 

and if {bo,uo) then (6 ^,#,Jq) converges weakly to {b,u,b) where {b,u) stands for 

the unique solution of 

dtb + u ■ Vb + ki ( 6 )div u = 0, 

dtu -t- u ■ Vtt — k 2 {b)Au + {ks^b) + n~^ki{b))Vb = 0. 

supplemented with initial data {bo,uo). 


(6.9) 


Proof: Let jf = jf + From (16.81) . we have 




-1 


n_i -I- u||V6^||^’^ n + n_i YCtJU. 

) lHr+AL) ) 


n_-I n 

SBI, +BI,) 


Hence, given p6.1l) . we deduce that 

( 6 . 10 ) jf = 0{£) in L^( 

Using the equation of jf, this gives 

(6.11) Vjg + £(1 -|- Cs)ji —>-0 in the sense of distributions. 

As in the non-equilibrium case, (16.8p implies that the families ( 6 ^) and (u^) are bounded in 
L°°(M+;H|i nH|JnLi(M+;.B|i +H|^) and L°°(M+; .624 ), respectively. 

Hence {dtb^) is bounded in L^(R+;ii 2 ^^ ) and we can thus deduce that there exists b in 

L°°{M.+ ; BfCiBfi) (b L^(M+; 112 ^^ + Bf^) and a sequence {£k)k&N going to 0 so that, for 

all (/) € 5 and all a G ( 0 , 1 ] 

( 6 . 12 ) cjb^'^^cjb in L°°(M+;H 2 y“). 

For {vS), we still have the weak convergence result given by (j5.37p . which suffices to pass to 
the limit in the mass equation. 
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Next, we observe that (16.8p implies that (jq) is bounded in B 21 +-^ 2^1 )■ Hence, 

there exists a sequence {ek)k&N going to 0 so that jg'' ^ jo in the sense of distributions. 
Moreover, we have 

jo -b^ = -C~^n~^diYjl - eCr^dijl- 

Remembering (|6.10p and (j 6 .ip . we see that the first term of the r.h.s. is 0(e) in a suitable 
space. The second one also tends to 0 in the sense of distributions as —>• 0. Hence 

(6.13) jo = b. 

In order to pass to the limit in the velocity equation, we proceed as in the non-equilibrium 
case. First we use (|5.37p and next, the fact that 

dt(ir + - k4mjA = -u^ ■ V# + k2mAiF - kAbAVW + -k'M)dtb^A - - kAbAVjL 

\ n J n n 

Taking advantage of (16.8p and of product laws in Besov spaces, we readily obtain that the 

four first terms of the r.h.s. are bounded in (M_|_; ) (or only in B 2 ) if 

n = 2). To handle the last term, we observe that according to (j 6 . 8 [l . Vjo is bounded in 

Tl r\ , O . ^ O 

B 21 ) hence the term A: 4 ( 6 ^)VjQ is bounded in ) (or B,^^ ) 

if n = 2 ). 

As in the non-equilibrium case, it is now easy to conclude that there exists some v in 
■ -—1 

L°°(M+; B 21 ) so that for all 0 in 5 and a G ( 0 , 1 ), we have 

h 

( F -* \ • IL — l—rv 

+ in L£,(M+;R|i ). 

Of course, combining with (I6.10p . this implies that v = u, and {b, u) thus satisfies the second 
line of (|6.9p . 

Finally, that the whole family ( 6 ^, jg, jf) converges to ( 6 , u, 6 , 0) stems from the fact 
that the solution to p6.9|) is unique (note that it is just the standard barotropic Navier-Stokes 
equations with a modified but still stable pressure law). □ 

7. The Poisson diffusion regime 
This section is devoted to the study of the asymptotics regime where 
(7.1) e <C T < and C^Cg ~ 1. 

According to the formal computations of Section [21 we expect the solutions of (|1.6p to tend 
to those of the Navier-Stokes-Poisson sytem (12.91) . 

The general scheme of the proof that we here propose is the same as in the study of the 
other asymptotics: we first perform a fine analysis of the linearized equations so as to check 
the long-time stability and exhibit the quantities that are likely to be bounded uniformly 
when e ^ 0, then tackle the proof of the global existence. We rapidly justify that the limit 
system is globally well-posed in a functional framework that is consistent with the analysis 
we used for m, and eventually take advantage of compactness arguments so as to prove 
the expected convergence result. As in the other regimes, the fact that the limit system has 
a unique solution will guarantee that the whole family of solutions to (II. 6 h converges to the 
solution to ( 12 .9p . 
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7.1. Linear analysis of (11.61) in the Poisson regime (17.111 . We here gather the estimates 
for (14.711 that have been obtained in Section 01 Recall that £ := vC. 


Small frequencies. Using (I4.28h . (j4.33|l . (I4.34h and the fact that \{b,d,jo,ji)\ « |(b, c), joJi)! 
and that the last term in the original definition of jo in (14.2211 has a negligible contribution 
with respect to ji, we get 

(7.2) \(b,d,^jo,pjo,ji){t)\ + f \(b,d,jo,ji)\dT + f \}o\dT + -[ |pCo| 

L JO Jo ^ Jo 

H-/ \h\dT < C\{b,d,^jo,pjo,ji){0)\ for all 0 < jO < Cl, 

with To := lo-V^h- ^^pd, Co := To - and ji := Ti - ;^To + 

nC‘^M 


Middle frequencies. Combining (j4.54p and the definition of Ci versus that of ji, we get 

(7.3) \{pb,d,jo,ji){t)\+p [ \b\dT + p‘^ [ \d\dT + p'^-f \jo\dT 

Jo Jo £ 40 

CM - _- 

H-/ |ji| dr < |(p5,d, jo,ji)(0)| for Ci < p < c£A4. 

£ 4o 

Large frequencies. Finally, using ()4.70p . we have 

__ ft ^ ft ^ CM ^ ^ 

(7.4) \{pb,d,jo,ji)it)\+p"^ \d\dT + p |6|drH-/ |(jo,4i)|dr 

4o 4o £ Jo 

<C\{pb,d,jo,ji){0)\ for p > cCJLi. 

Therefore, localizing (14.4p (with nonzero source terms / and g) according to Littlewood- 
Paley operator A^, using (14.Sp . following the computations leading to the above three in¬ 
equalities and using Fourier-Plancherel theorem, we end up with the following inequality for 
all s E M 


(7.5) ii(£,ji)(t)iis,+ii6(t)iT.i +\\bm^i, + i\\jom% +ibo(t)ir-(+i + iijo(oii 


ibi 




h,l 

m 


+ 


/ ||tt|| os+2 dr + 

Jo 2,1 


+ 


CM 


2,1 
Zm 


^Jo,ji)|||!+2 dr + ^ 


/ \\h\fZ" 

Jo 2.1 


dr + 


£ 


m,l,£A 4 I 

ns+2 I 


'0 


>0 


f IIColl“+. 

Jo 2,1 


l/l,l I I 

I 0.+1 dr + 


dr + / 

J{ 

CM 


kWpl dr 
-^2,1 


e 


[ IKioJOII^f'^dr 

Jo 2,1 
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< iiuJi)(o)iIb. , + iii-mic; + iii.(o)ii*..‘„ + iiuown'.; + niomir;*. 

2.1 ^ 2,1 ^ 2,1 C 222,1 122 1 


with jo ■■= jo - Vnb - y/n J div£, 


Co ■= jo - y/nild- 


1 


nC^M 


+ I bo ( 0)11 + 


-1 


£,i 


B. 


+ 


2,1 


iSs+l 

-^2,1 


+ lldb|Jdr, 


A b and ji := ji + 


Vjo - —V6. 


y/nCM 


CCsM 
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As in the previous sections, owing to the convection term in the equation for b, the above 
inequality does not allow to prove the global existence for (II.61) . and one has to consider the 
paralinearized system (I5.10I1 . Adapting the proof of Proposition 15.11 we get the following 


Proposition 7.1. If the coefficients £, Cg and e fulfill (17.ip then for any smooth solution 
{b,u, jo, ji) to (|5.10l] . one has the following inequality 

\\{b,u,jo,jl){t)\\Y.^t^^ < c(^\{b,u,jo,ji){t))\\x. + \\Vv\\L°°\\{h,u,jo,ji)\\x. dr 


+ 


||(VF,G)||H_,dr+ f\\{F 
^ 2,1 -^0 


-T^-Vb,G-T^- Vu )\\^'\_i dr 


^2,1 


with 




B 2 
^2,1 


R 2 
^2,1 


It- I ^ L/U _ 

B^ Lv B^ 
^ 2,1 -^ 2,1 


B 2 
^2,1 


B 2 
^2,1 


and 


ji)(i)||y;^(i) := sup \\{b,U,jo,ji){T)\\x^ 


0<T<t 


[ (II^L^+1 + Il(ftji)ll%+1 + +J [ \\Co\\%'dr 

Jo ^2,1 Rj 1 ^2 1 ^ Jo ^2 1 


+p 


-1 


ioll . dr H- 


'0 ^2,1 


nfi^,dr + 


'0 ^2,1 


/ \\bt\ 

'o 


dr + 


CM 


ll(jo,ji)||'*’S'^^ dr. 


Bo 


Above, we set 


io := io-ft-^divu, Co:=jo-(Id- 


6 and 


7.2. Uniform global well-posedness in the Poisson regime. In this paragraph, we 
sketch the proof of the following global existence result. 


Theorem 7.1. There exists a positive constant rj depending only on p-jv, n and on the 
pressure law such that if e G (0,1) and if the coefficients C and Cg fulfill ()7.ip then any data 
(^ 0 )^ 0 )io, 0 ) Ji,o) satisfying 


C^-6) ll(^05'“0)i?0,05il,o)llA|' ^ 

generates a unique global solution {b^,u‘^, Jq, Ji) Yf' to System (11.611 . 

Furthermore, we have 


(7.7) 


Wib’^G ogDWyj' < C\\{bl,il^, jo 0, ji o) 


lAr • 


Proof. Assuming with no loss of generality that p = 1, the proof relies on Proposition 17.11 
with, dropping the indices e for better readability, v = u 

CM 

F :=—Txh ■ u — ki{b)div u and G :=-T^^-u + k 2 {h)Au — ko{h)Vh-\ - ki{b)ji. 

Let us just explain how to handle the last term, as it cannot be bounded exactly as in the 
proof of Theorems 15.11 or 16.11 due to the difference between the spaces Yf and Yf. We use 
the fact that 


CMh = CMh + CjWb - Vjo, 
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and thus 

CMh{b)ji = + C-^ki{b)W^^^ - 

It is clear that 


wkmji 


-th,CM I 




< 


L°° 


(4^1) 


-^h,CM I 

h \ 


L^(B? 


< 


{CM) \\ib,u,jo,ji)\\ 




2,1 


that the second term in the r.h.s. may be bounded in the same way, and that the third one 
can be bounded as the pressure term k^{b)Vb. For the last term, one just has to observe that 
the definition of || • ||yi guarantees that 


(7.8) 
and that 


l,CM 


llVjoll ^ 
11^4(6)11 


< 


\\{b, u,joJi)\\ 




n < 


The rest of the proof is standard, and thus left to the reader. 


□ 


7.3. Study of the limit system. We introduce the following norms 


\\{b,u,jo)y, := 


,e,u- 

n 

SoT 




I 

n 

b7, 


+ llf^l 


B^~ 

-^2,1 


+ 


n 

b7. 


+ V'- 


bI 




and 


\\{b,u,jo)\\y, := sup||(6,n,jo)(i)||;e- 

t>o 


+U 


/ (ll('>.io)ll 

J M_l 




^2,1 


+ 1 


+ \\u\ 


B^ 

^2,1 


+1 


) (it + f 

Jr. 


h,u 

n 

B^ 

^2,1 




n 

B^ 

^2,1 


+2) dt- 


Theorem 7.2. Let the data (f>o, f^o, Jo,o) satisfy for a small enough constant c > 0 


(7.9) 



+ ^^ll^o|| 


h,u ^ 

n 

Bh 


^2,1 


and the compatibility condition 

Jo,o-^Jo,o = oo- 

nm 

Then System (1^ admits a unique global solution (5, ft, jo) in the space satisfying in 
addition for a large enough constant C independent of v 


(7.10) 


jo)||y- < + u||6o||^’? 

D 2 D • 

tSn 1 XJo 



)• 


Proof. We just sketch the proof as it is very similar to the standard one for the barotropic 
Navier-Stokes equations. As usual, it suffices to treat the case u = 1. 

The first step is to analyse the linearized system 


' dtb + divft = /, 

< dtu - Au + Vb+ 1 Vjo = g, 

- (Id - = ft- 


(7.11) 
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To this end, we set d = (—A) ^/^divu and observe that in the Fourier space, {b,d) fulfills 
the following ODE f = g = 0 


dtb + pd = 0, 

dtd + p^d — pupb = 0 


with p := 1^1 and Op := 1 + 


Of course, jo may be computed from b by the relation 


^ nm y 

Jo = - b. 

p^ + nm 


Introducing the following Lyapunov and diffusion functionals 


£p = 2ap|6p + 2|(ip + |p6|^ - 2Re (pftd) and Til = p‘^{ap\b\‘^ + \d\‘^), 


we see that 


1 d o o 

2^4+«! = »• 


Because we have < cq for some cq independent of p, one can thus conclude exactly as 
in the standard barotropic case that for all t > 0 and p > 0 

\(b, pb, d){t)\ + mm{l, p) f \pb\dT + p'^f \d\ dr < \(b, pb,d){0)\. 

Jo Jo 

Back to (17.111) . one may combine Fourier-Plancherel theorem and Duhamel formula to get 
the following estimate for all s G M 

II(6, V6,tt)(t)ll^s + Ib'o(f)llw [ (ll’^ll_Bf+^ “ 1 “ IK^jJo)ll^s+ 2 ) 

2,1 2,1 2,1 Jq 2,1 .02,1 

rt 


+ 


[ (ll?^ll 5 ’.Vi + l|jo||^’l+3)^iT<cf||(6o,V6o,«o)|lB| + / liav/,; 

Jo " 2.1 " 2,1 V Jo 


dr . 


2,1 


However, because of the convection term in the equation for b, this does not allow to prove 
estimates for the nonlinear system (12.9h . Therefore, mimicking the standard approach for 
the compressible Navier-Stokes equation we ‘paralinearize’ the system and get the following 
result 

Proposition 7.2. The solutions to the following paralinearized system 

dtb + Tjj ■Vu + div u = f, 

' dtU + n-Vu-Au + Vb+j^Vjo - g, 

fulfill the following a priori estimate 

||(6,V6,u)(t)||^. + IIjoWIIb. + [ {\\u\\B^+fi + \\{b,jo)\ffil+ 2 ) dr 

2,1 2,1 2,1 Jq 2,1 - 02^1 

+ All^ll^’.V + IIJo||^i3)dr<cf||(5o,V6o,no)|ls. + ril(/,V/,5)||^. dr 

JO 2,1 -“2,1 V 2,1 Jq 2.1 
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Now, in order to estimate the solutions of the nonlinear system (12.91) . it suffices to apply 
the above proposition with v = u 

/ = —• V6 — A:i(6)div ft and g =—T!^-Vu + k 2 {b)Au — k:i{h)\/b — n~^ki{b)VJ q. 

All the terms but the last one of g are already present in the barotropic Navier-Stokes 
equations, and may be bounded quadratically in terms of ft, Jo)||j)i- Now, we have 


11^4(6) VjolL: 


- 1 < 


^2,1 


^2,1 


hence 


||A:4(6)Vjo|L,.t-t 


< 


^,1 


,(lboir;. 


+ 


'L\{BIA) - 

and one can thus conclude that whenever the solution {b,u,jo) exists we have 

which allows to get (|7.10p if ()7.9I) is fulfilled with a small enough c. 


□ 


7.4. Weak convergence. Here we justify weak convergence to (12.91) when assumption (17.11) 
is fulhlled and, in addition 

(7.12) ^ m € (0, +oo). 

Theorem 7.3. Let the family of data {bQ,UQ, jo di o)o<e<i satisfy (|7.6p . Then the global 
solution (6^, ft^, Joiii) ^e' given by Theorem \7.1\ satisfies 

(7.13) j! = Oi£) m +i?|i), 

and, up to extraction, {¥,hT,jff) converges weakly to some solution {b,u,jo) in of System 
(EJI) when e goes to 0. 

If in addition 

(7.14) {bo,ifo,jo,o) {bo,uo,jofl) with - u^Ajo,o + R"i(jo,o - bo) = 0, 

then the whole family {b^,vF,jQ) converges to the unique solution {b,u,jo) corresponding to 
the initial data {bo,uo, jo,o), given by Theorem \ 7. 4 

Proof: Let us first prove ()7.13p . From (17.71) . we already know that (jf)^’'^'^ and 
are 0{sC) in Bfi )• Now, we have 


jf = i! 


-vjg + 


-VbS 


CCs ^ CCsM 

It is easy to see that the last term is 0{C^) in Bf^ + Bf^ ), and that, according 

n 

to (I7.8p and C^Cg ~ 1, the last but one term is 0{C) in Bfi), which completes the 

proof of (j7.13l) . 

Next, let us turn our attention to the convergence of Jq. First, (17.71) and the definition of 
II • lly,. ensure that (jo)^’"^'^ is bounded in, say, ). Next, using the bound for 

the middle frequencies of jo and for the low frequencies of ^o, we discover that (jo)^’"^'^ is 

n 

bounded in L^(M+; i? 2 ^^). Hence, up to an omitted extraction 


Jo 


jo weak * in L^(R+;ij 2 ^^ +Bf-^^). 


(7.15) 
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Now, taking the divergence of the equation of jf, then using the equation of jg gives 

(7.16) = -C{1 + Cs){nC{y' - jg) - endtjo) - edtdivjf. 

Given (|7.13p . one can assert that the last term tends to 0 in the sense of distributions. We 
also know that, up to an omitted extraction, jg —>■ jg in the sense of distributions, hence 
given that £(l + £s)e^0, the term with dtj^ also tends to 0. Finally, exactly as in the 

, n_2 . IL 

cases treated before, (6^) is bounded in L°°(M+;i72^^ H hence weakly converges to 

.71 _^ .72 

some b G (R+; 7^2^ n B^i). As (17.120 has been assumed, passing to the limit in (|7.16p 
gives 

i^^Ajg = -nm(b-jo). 

Passing to the limit in the equation of b goes along the lines of the non-equilibrium case 
we notice that (dtb^) is bounded in B21 ) and we thus have, up to an omitted 

extraction 

72 _ 

(7.17) (jb^^cjb in L°°(R+;5|i ) for all a G (0,1). 

As (j7.7p also implies that (n'^) is bounded in L°°(R+; B21 )nA^(R+; B21 ), we have uF ^ u 
weakly * in that space, which is enough to justify the first equation of (|2.9p . 

In order to pass to the limit in the velocity equation, we use again the fact that 

dt(if + - hib^^) • V# + k2{W)Av^ - + -k'M)dtb^ Jf - - ki{¥) VjL 

\ n j n n 

As in the other asymptotic regimes, the first four terms of the r.h.s. are bounded in 
L^(R+;i72^^ ) (or in B^^o ) if u = 2 ). To handle the last term, we observe that 

according to (j7.7p and (I7.15p . (Vjg) is bounded in L^(R+;i72^^ + B21 ). Because {¥) is 

. ^_1 . 21 .21 _9 

bounded in B21 FiB^i), this implies that k4 ^{¥)VJq is bounded in L^(R+;i?2^j^ ) 

(or L^(R+; B20Q ) if ti = 2 ), and thus ^ ki{¥) jf) is bounded in the same space. 

As in the already studied cases, we conclude that there exists some u in L°°(R+;i?2^i ) 
so that for all in 5 and a G ( 0 , 1 ), we have 

^ in T^c(I^+; ^2^^”“)- 

Finally, in the case where ( 17 . 141 ) is fulfilled, the limit system (| 2 . 9 p supplemented with initial 
data (60) jo,o) possesses a unique solution (6, ft, jg) given by Theorem I 7 . 2 P and the whole 
family {¥,u^,jo) thus converges to {b,u,jo)- Gl 

Appendix A. Estimates for a toy linear differential equation 

The appendix is devoted to the proof of decay estimates for the solutions to systems of 
ODEs of the form 

(E) dtU + AoU + piAi + Bi)U + p‘^A 2 U = 0 , 

where p is a nonnegative parameter, and Ag, Ai, Bi and A2 are given N x N matrices. 
We have in mind System (14.7p in which case, after suitable change of unknowns (see (IA.4p i. 
Ag is a degenerate nonnegative diagonal matrix, A2 has nonnegative eigenvalues and Ai is 
skewsymmetric up to some positive diagonal symmetrizer. 
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A.l. A general approach. The basic idea is to set V := (I + pP)U where P is a suitable 
matrix, so as to eliminate the bad first order term pBiU. Now, whenever {I+pP) is invertible, 
the equation for V reads 

dtV + AoF + p{Ai + Pi + [P, Ao])V + p^Ao, P]P + [P, Ai] + [P, Pi] + Aa)^ 

+pHI + pP) ((Ai + Pi)p2 - AoP^ - A 2 P) (I + pP)-^V = 0. 
Therefore, if one can find some matrix P so that 
(A.l) [Ao,P] = Pi, 

then we have 

(A. 2 ) dtV + AoF + pAiV + (A 2 + PPi + [P, Ai]) V = p\l + pP) A 3 (/ + pP)"V, 

where A 3 := (PAq — Ai)P^ + A 2 P 

The gain is clear as the matrix Pi now appears at order 2 instead of order 1. Hence the 
system for V is more likely to be tractable for small enough p as we shall see below. 


A.2. Application to the linearized system for barotropic radiative flows. The sys¬ 
tem we are interested in reads 

l^\ 

0 

0 ’ 

VO/ 

To bound the solutions of (IA.3P for small enough p (under some stability condition that 
we will discover below), we propose two different approaches, the first one being appropriate 
to handle the case where (5 and 7 are of the same order of magnitude, and the second one, 
to the case where / 3 / 7 ^ 1 or 7 // 3 <Sl (of course only 7 > /3 is relevant as far as (14.71) is 
concerned). 




( ^ 

\ 


/ 

0 

P 

0 

0 

V 


( ^ 

\ 

(A.3) 

d 

d 


+ 


-p 

P^ 

0 




d 



dt 

jo 




-r] 

0 

/3 

ap 



jo 




\ ji 



V 

0 

0 

—ap 

7 

) 


\ ji 

/ 


where all the coefficients of the matrix are positive. 


A.2.1. First approach. Making the change of unknown 

/ 1 0 0 
0 


(A.4) 


U := 


\ 0 


0 

1 

0 

0 


0 ^ 

7 

1 0 
0 1 


\ 

/ a \ 


d 


h 

/ 

\ ji / 


and setting a := a + we see that U satisfies a system of type (P) with 


/ 

0 

0 

0 

0 \ 




/ 

0 


1 

0 

0 

\ 


0 

0 

0 

0 



Ai ■= 


1 - 

a<;'n 

^7 

0 

0 

0 



0 

0 

/3 

0 





0 


0 

0 

a 


V 

0 

0 

0 

7 




V 

0 


0 

—a 

0 

/ 


/ 

0 

0 

0 

£ 

7 

\ 




( ° 

0 

0 

0 

\ 



0 

0 


0 



A2 


0 

1 

0 



_ 


0 


7 

0 


and 

• — 

7 




V 

13 

0 




0 

0 

0 

0 



V 

ar) 

T 

0 

0 

0 

/ 




V 0 

0 

0 

0 

/ 
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Note that the above matrices may be written in block form as follows 


= 


0 B\ 
0 


Bl 


An = 


0 0 
0 A 


^1 = 


_M} 0 


0 Af 


P = 


pll pl2 

p21 p22 


Computing the commutator 

(A.5) [Ao,P] = 


0 -P^'^A 

AP21 [A, P22] 


we see is satisfied if 

pll := 0, p22 ;= 0, := -Pi A-\ P^i := A~^P?. 

In other words 

/ 0 0 0 

. s 0 ( 

(A.6) P= Q 

0 

V 0 0 0 

which, remembering (IA.4() . corresponds to the following change of unknowns 


\ 


^ 0 
0 0 


/ 


(A.7) 


V = 


/ b \ 
h 
jo 

vTi y 


/ 1 


— R 

P 


0 0 
1 


7 


rP 


\ 


-WP 


-^P 0 


/3^p 7 

1 -^P 

P 


\ -T,P 


0 


1 



/ a \ 


d 


Jo 


\ ji / 


Note that the determinant of the matrix (I + pP) is 




1 2 


and is thus of order 1 whenever p satisfies the smallness condition 


(A.8) 


o2 < 




min(/3^,7^). 


In order to go further in the estimates of V, we compute 


PPi = 


-piA-iPi2 

0 


/ _fA2. 


A-iPi2p| 


0 


\ 


0 

0 


0 

0 


0 

0 

0 


0 \ 
0 
0 

aqr) 

WA } 


iPAI-f 0 -BlA-^Aj + AlBlA-^ 

[P, AiJ - _ ^2^-152 0 


/ 0 

0 

gg’? V (-1 

Pi ^ 1 p-f' 

V 0 


0 

0 

0 

— RR(L-^l\ 
P ^P^i> 


_“£('!+If 

7 ^ /3 ^ 7 ^ 

0 

0 

0 


and 


gg<? \ 

Pi ' 1 ^ V / 


Pi 


0 

0 
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Finally, ^3 := {PAq — Ai)P'^ + A 2 P reads 
(A.9) 


A3 = 


a?r/ 


( ^ ^ 

i — 1 h 4-^1 n i — 

7 ^ \ /Sy ) 77 

0 


V 


0 


\ 

^ 0 
V 

0 0 ^ 

0 0 


Therefore, resuming to (IA.2h . we conclude that 

jV + AqV + pAiV + p^{PBi +A2)V = p‘^[Ai,P]V + O(p^). 

Of course, the remainder term 0{p^) strongly depends on the coefficients of the system. 
We shall see below that the structure of [Ai,P] will enable us to treat p^[Ai,P] and the 
nondiagonal term of A 2 as small error terms as well. 


Let us focus on the system satisfied by (b, D) for a while. We have 




+ P 


0 

PI 


+ P" 


= P 


“iii ii 

7 ^ /3 ^ 7 / 


WF 

0 

0 


0 

1 _ 

Wi 


_ ggy _ ? n 4- 

'Y Pj B'Y 


Pi 


h 

b 

jo 


+ o{p^). 


For small enough p, optimal estimates may be proved by taking advantage of the results of 
Appendix [HI Indeed, denoting by Fp the r.h.s. of (jA.lOjl . we see from (IB.6p that if we set 


Ui-= [I 




then, under the following necessary and sufficient stability condition 
(A.ll) 

we have (see (IB.4p and (jB 
(A.12) U. 

whenever 




p~|(b,b)| and jU^^ + ’^p‘^U^^<Up\Fp\, 


(A.13) 


P < 


1 + ^ 

Pi 


1 _L ^(-1 _ 1 ) 

^ + /37 W P) 

So finally, we get for some appropriate constant C = (7(0;,/3,7, ?, ?/) 

\{b,d){t)\+9p'^ |(b,0)|dr < (7^|(b,'D)(0)| |(jo, ji)| dr + |(b,b, jo,ii)| dr^ , 


which, if p <C f', may be simplified into 


(A.14) 


|(b,b)(t)| +i>p' 


Jo 


dr < Ci |(b,c))(0)| +p' 


ri(jo,Ti)| 

Jo 


dr 
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The modified radiative modes jo and ji fulfill 


I t hK- 0 


0 


I » 7 + lfp^ 


h 


jo 

h 


= P 


aar) ^ M _L 
0 


0 

^(L _i_ i'l 
/3 1/3 7I 


+ o{p^). 


Therefore we easily get 


14 (|T„|=+“fcid + (^] iToi= iti^ 


2 dt 


a 




a 


< C(p^|(b,h)| +/9^|(b,h,jo,ii)|)- 

Then, integrating and assuming that p <C 1 yields 

(A.16) (|0o,ji)(t)| +min(/ 3 , 7 ) / |0o,h)| dr < C'flOo, ji)(0)| + / |(b,b)|dr 

Jo V do y 


Combining with (jA.14p , we can conclude that there exists some positive constants po and C 
depending only on (a, /3, 7 , p) so that for all 


(A.17) 
we have 


0 < p < min(l, D) po, 


(A.18) |(b,b)(t)|+p|0o,ji)(t)|+i>p^ / |(b,b)|(ir+ pmin(/ 3 , 7 ) / |(jo,ji)|dT 

Jo Jo 

< C(|(b,b)( 0 )| +i>|Oo,ji)(0)|). 


A.2.2. Second approach. In the case where /3 and 7 are not of the same order of magnitude, 
Inequality ()A.18p is not fully satisfactory, first because we would like to have a control on 
(3 ffl |jo| dr and 7 |ji| dr rather than just on min(/ 3 , 7 ) |0oJi)| dr and, second, because 

the range for which (IA.18P holds true tends to shrink to 0 if /3 <C 7 or 7 <C /3. 

In this paragraph, we propose another approach to handle (IA.3I) in the case /3 7 ^ 7 , still 
based on rewriting the system in the form ()A.2p . but with a different definition of Ai and 
Bi (Aq and A 2 being unchanged). More precisely, we now set 


/ 


Ai : = 


1 0 0 \ 


-l-f2 0 0 0 


V 


0 

0 


000 
000 / 


and Bi := 


( 0 
0 

V 


0 


0 


0-^0 
0 0 a 

^ 0 -a 0 y 


Then writing the matrices coming into play in block form, we see according to (IA.5I) . that a 
possible choice for P is 


pii := 0, 




1 


P-1 


0 d 
a 0 


with d := a + —— 

Pi 
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With this new definition of P, we have 

7- ° 

0 

aari 


I 


PBi = 


CX.OLT] 

0 


a<;ri 

Wi 

0 




0.7) 


0(^7) I OO 

Wl 7-/3 


0 

oo<; 

Pi 

0 


^ _i_ 

/ 


and [P, ^i] = 


_ OO 

' P-i / 


0 

0 


■(1+s) 


Therefore setting 
(A.19) V = 


( b \ 
0 
io 

/ 


1 

Wi^ 
— R 

/3 

. 07) 

\ lii-P)^ 


0 

1 

-jiP 

0 


0 

1 

■^P 


0 

0 

0 

_ 07) 

Pi 


-RS. 

Pi 

0 

0 

0 




0 \ 
( 1 +w) 


0 

0 


/ 


\ 


f-i_ 

\/3—7 p^"j) P 

1 


/ 



/ “ \ 


d 


jo 


\ ji / 


it is clear that working with {a,d,jo,ji) or (b, 0 ,jO)ii) is equivalent whenever 

P<Ch-/3\, 

for some positive constant C depending continuously on the coefficients of the system. 
Putting together the previous computations, we see that V fulfills 

/ 0 0 0 \ 

0 0 0 


A. 

dt 


F + 


V 


0 

0 


0 

0 


0 


0 


V 


+p 


= p 


0 

0 


l+m+J^)p^ 

0 1 0 0 \ 

V 

0 0 0 
0 0 0 / 

(i+i) 


-1-i^ 0 0 0 

Pi 


007) 


V 


Pip-l) 


a_ V (-] iRSR\ 

-l) W Pi ! 


0 
0 

0 

0 

0 

07) 

lip-l) 


“i _ _ 

1 Kp^ 1 

0 


actj ? ll I 


Pi 


0 

0 


0 

0 


/ 


+p\l + pP)A3iI + pP)-% 

with A 3 := {PAq — Ai)P‘^ + A 2 P satisfying IA 3 I <C{l + \^ — /3\~^). 

Next, arguing exactly as to handle (lA.lOp . we discover that under the stability condition 
(|A.lip and for p satisfying (jA.lOp (and of course also p < c|/3 — 7 ^), we have 


(A.20) \{b,d){t)\+Dp'^ |(b,'D)|fir<C' |(b,0)(O)| 


+ p- 


Jo 


dr+ p^[l + 


It -/^l' 


|(b,0,jo,ji)|dr 
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Now, in contrast with the first method, we can bound jo and ji independently from one 
another from the equation satisfied by jo, we readily get 


(A.21) |jo(t)| + ( /3 + liol dr < |jo(0)| + ^ 1^1 dr 


+ Cp^ 1 + 


1 


|7-/3|3 


|(b, 0 ,jo,ji)|dT, 


while ji satisfies 

(A.22) + + + 


aa 


Cp^ /■* r, j 

+ -T / lt)| dr 


+ Cp^{ 1 + |^_P| 3 ) ^ l(^5,io,ji)|dT. 


Putting inequalities (IA.2nh . (IA.2ip and (IA.22P together, it is now easy to conclude that 

rt 


(It 


(A.23) |(b,h)(t)|+z>|7-/3||(jo,ji)(t)|/ |(b,b)| 

Jo 

+ p| 7-/3|^/3^ |jo|dr + 7^ |ji|dp < C'(|(b,b)(0)| + i>|7 -/3||(jo,ji)(0)|), 
if, for some small enough constant c depending continuously on a, j3, 7 and ?, we have 
(A.24) p<cminfl, p, | 7 -/ 3 | 2 , i>| 7 -/ 3|3 


(B.l) 


Appendix B. Optimal decay estimates for a toy system 

For the reader convenience, we here recall some results that have been obtained in our 
recent work [ 8 ] for the following linear system of ordinary differential equations 

dtX + apY - bp^X = A, 
dtV - cpX + dp'^Y = B. 

Above, p stands for a given nonnegative small parameter and a, b, c and d are four real 
numbers satisfying the stability condition 

(B.2) a > 0, c > 0 and d — 6 > 0. 

Routine computations show that the following Lyapunov functional := c\X\^ + a\Y\‘^ — 
p{d + 5)Re (XY) satisfies the relation 

(B.3) 147 + + (fl^)7Re(AF) 


= Re {cAX + aBY - p{b+d){BX+AY)). 


Now, observe that whenever p < we have 


62-^2 


p^Re (XY) 




pc|X|2 + a|y|2) < £2 < + a|y|2). 


and 

(B.4) 
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which leads if A = i? = 0 to 

(B.5) 

and thus 

(B.6) Cp{t) < 

Combining with (jB.4p and Duhamel’s formula, we deduce that for general source terms A 
and B we have 

(B.7) Vc|X(t)|2 + a|y(t)|2 < V3e-(^)^'‘(^v'c|X(0)|2 + a|y(0)|2 

+ f y^c|^| 2 + a|-B|2 dr 

Jo 
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